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ABSTRACT.  A  technique  for  determining  the  mass  of  large 
rocket  motors  during  burning,  and  an.  experimental  system 
developed  for  use  with  the  NOTS-designed  three- component 
static-test  stand  at  the  Skytop  facility,  are  described. 
The  technique  involves  mounting  the  rocket  motor  on 
springs  and  continuously  exciting  the  spring-mass  sys¬ 
tem  at  its  natural  frequency. 

The  theoretical  design  of  the  control  system  is  dis¬ 
cussed  and  theoretical  system- performance  characteristics 
are  compared  with  those  of  the  experimental  system.  Re¬ 
sults  of  an  analog  simulation  study  conducted  in  two  de¬ 
grees  of  freedom  are  included  in  the  appendix. 

Analysis  of  data  from  an  evaluation  test  involving  the 
static  firing  of  a  large  rocket  motor,  indicates  that  a 
2$>  accuracy  in  mass  measurement  is  attainable  with  the 
experimental  system.  Criteria  for  the  design  of  an  oper¬ 
ational  system  that  will  produce  an  accuracy  of  1$  are 
presented. 
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FOREWORD 

In  the  static  testing  of  rocket  motors,  the  prime  requirement  is 
to  he  able  to  determine  accurately  the  force  of  the  developed  thrust. 

With  the  advent  of  large  rocket  motors  it  became  necessary  to  test  them 
with  correspondingly  large  static-test  stands  having  very  low  natural 
frequencies.  However,  because  of  the  low-frequency  characteristics  of 
the  test  stand,  oscillations  occurred  during  irregular  burning,  thrust 
buildup,  and  termination  periods  and  these  extraneous  oscillations  were 
recorded  as  part  of  the  thrust  load-cell  data.  To  separate  oscillation- 
caused  data  from  the  thrust  data,  a  data  combination  system  was  developed., 
To  use  the  system  it  is  necessary  to  know  the  mass  of  the  oscillating 
system  throughout  rocket -motor  burning. 

This  report,  issued  at  the  working  level,  describes  a  technique 
for  determining  the  mass  of  a  system  continuously  during  rocket-motor 
burning.  It  also  includes  a  description  of  the  theoretical  development 
of  an  experimental  mass -measuring  system,  the  design  and  evaluation  of  a 
prototype  system,  results  of  an  analog  simulation  study,  and  information 
on  full-scale  testing. 

This  project  was  started  in  the  spring  of  1959  under  Task  Assignment 
SP  JlkOl -7  and  was  completed  with  an  evaluation  firing  in  December  1962 
under  Task  Assignment  SP  71402-8.  The  report  has  been  reviewed  for  tech¬ 
nical  accuracy  by  Benjamin  Glatt,  R.  W.  Murphy,  D.  P.  Ankeney,  and  C.  E. 
Woods . 
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INTRODUCTION 


In  the  static  testing  of  large  rocket  motors,  early  test  stands 
often  employed  a  standard  load  cell  to  solve  the  following  simplified 
thrust  equation: 

T  =  KX  (1) 

where 

T  =  thrust 

K  =  spring  constant  of  load  cell 

X  =  spring  displacement 

It  soon  "became  apparent,  however,  that  the  necessarily  low  resonant- 
translation  frequency  generally  associated  with  large-scale  testing  re¬ 
sulted  in  the  occurrence  of  undesirable  oscillations  or  'ringing*  during 
thrust  buildup  and  thrust  termination.  This  is  the  result  of  the  thrust 
buildup  and  termination  appearing  as  approximate  step  inputs  to  an  under - 
damped  stand  load-cell  system.  To  account  for  these  oscillations,  the 
basic  single  degree  of  freedom  equation  can  be  used. 

Consider  the  rocket  motor  and  stand  illustrated  in  Fig.  1 


FIG.  1.  Rocket  Motor  and  Thrust  Stand  Schematic 
Illustrating  Basic  Thrust  Equation. 
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Where: 


T 

K 

C 

M(t) 


c.g. 


X 


developed  thrust,  lh 

spring  constant  of  load-cell,  lb/ft 


lumped  damping  constant, 


lb-sec 

ft 


mass  of  total  moving  structure  (varies  with  time), 


lb-sec 

ft 


center  of  gravity  location  (assumed  fixed  throughout  burn- 
ing) 


displacement  of  mass  center  of  gravity 


By  summing  the  forces  along  the  translational  axis  the  following  equa¬ 
tion  is  obtained: 


T  =  KX  +  CX  +  M(t)X 


(2) 


It  is  noted  that  this  basic  equation  is  a  differential  equation  of 
the  second  order.  As  such,  with  the  damping  constant  C  relatively  small 
(as  it  should  be  for  accurate  determination  of  thrust  buildup  and  termi¬ 
nation),  a  step  input  or  sudden  change  in  thrust  will  result  in  large 
overshoots  and  relatively  undamped  oscillations  in  X.  These  oscilla¬ 
tions  are  recorded  as  part  of  the  load-cell  data  and,  because  of  their 
low-frequency,  they  cannot  be  filtered  without  sacrificing  data  resolu¬ 
tion.  As  a  solution  to  this  problem  NOTS  has  considered  the  development 
of  a  data -combi nation  system  which  enables  the  separation  of  these  oscil¬ 
lations  from  the  load-cell  data.*  With  this  system,  data  from  an  accel¬ 
erometer  mounted  with  its  sensitive  axis  along  the  rocket -motor  line  of 
thrust  will  be  combined  with  the  standard  load -cell  data  and,  by  assuming 
a  linear  change  in  rocket  mass  during  burning  (as  determined  by  weighing 
the  rocket  motor  before  and  after  burning  or  as  calculated  from  internal 
ballistics  data)  the  basic  thrust  equation  can  be  solved.  The  acceler¬ 
ometer  data  is  integrated  to  obtain  X  and  by  estimating  C  the  second  term 
of  Eq.  (2)  can  also  be  included. 

While  in  many  routine  tests  the  assumption  of  a  linear  change  in 
mass  during  burning  might  be  useful  in  applying  accelerometer  data  to 
obtain  accurate  thrust  determinations,  it  is  not  adequate  for  use  with 


*U.  S.  Naval  Ordnance  Test  Station.  A  System  for  Correcting  for 
Spurious  Natural -Frequency  Ringing  of  Rocket  Static  Thrust  Stands,  "by 
J.  S.  Ward.  China  Lake,  Calif.,  NOTS,  1  September  i960.  (NAVWEPS  Re¬ 
port  7569,  NOTS  TP  2541) . 
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propulsion  systems  in  which  the  propellant  mass-burning  rate,  intention¬ 
ally  or  otherwise,  deviates  seriously  from  a  constant  value  during  burn¬ 
ing.  Furthermore,  it  is  often  desirable  to  know  the  mass-burning  rate 
to  the  same  degree  of  accuracy  as  that  of  the  thrust  at  all  times  during 
burning  in  order  to  evaluate  the  changing  effects  of  erosive  burning, 
combustion-chamber  free  volume,  operating  pressure,  charge  deformation, 
specific  thrust,  nozzle  erosion,  propellant  composition,  etc.  Likewise, 
it  is  very  important  to  have  continuous  mass  measurement  in  case  of  me¬ 
chanical  failure  in  the  propulsion  system  since  only  limited  specific 
thrust  and  specific  impulse  data  can  be  obtained  from  an  incomplete  test 
if  the  mechanical  condition  of  the  nozzle  and  the  amount  of  propellant 
burned  prior  to  the  failure  are  not  known. 

In  reviewing  this  problem  it  was  judged  that  a  reasonably  sophisti¬ 
cated  evaluation  of  test  records  would  provide  a  significant  increase  in 
information  if  a  system  could  be  devised  for  determining  rocket  mass  con¬ 
tinuously  during  a  firing.  With  the  high  cost  of  testing  large  motors, 
this  increased  information  would  create  substantial  savings  in  both  time 
and  money  during  developmental  programs. 

Two  methods  which  have  been  used  for  determining  the  mass  of  a  burn¬ 
ing  rocket  motor  are  referred  to  as  the  Z  Load  Cell  Method  and  internal 
ballistics  or  Pressure  Method. 

The  Z  Load  Cell  Method  involves  the  use  of  load  cells  in  the  thrust 
stand  legs  to  measure  the  vertical  force  during  burning.  If  the  stand 
and  thrust  misalignments  are  known  and  assumed  constant  throughout  burn¬ 
ing,  a  determination  of  mass  versus  time  can  be  made.  The  overall  accu¬ 
racy  in  mass  determination  is  estimated  as  however,  it  is  felt  that 
improved  stand  alignments  and  data  reduction  techniques  may  improve  this 
figure.  Continued  investigation  into  this  method  is  planned. 

To  use  the  Pressure  Method  of  determining  mass  requires  that  the 
rocket  motor  chamber  pressure  be  known  throughout  burning.  It  is  also 
necessary  to  know  the  total  amount  of  propellant  burned  and  the  throat 
area  of  the  nozzle.  Using  these  parameters,  calculations  of  mass  burn¬ 
ing  rates  are  made.  It  has  been  estimated  that,  due  to  the  difficulty 
in  measuring  these  parameters  and  the  linearizing  assumptions  made  in 
converting  the  information  to  mass-burning  rates,  an  accuracy  of  5$  is 
the  best  attainable  with  this  technique,  too. 

Using  this  as  a  general  guideline  in  determining  the  feasibility 
and  usefulness  of  a  particular  system,  the  following  technique  for  the 
continuous  measurement  of  mass  has  been  developed.  The  horizontal  NOTS 
Model  II  thrust  stand  used  in  Bay  I  at  Skytop  (Fig.  2)  was  selected  for 
system  development,  *  since  it  lends  itself  favorably  to  this  technique 
of  mass  measurement  and  since  considerable  knowledge  had  been  gained  of 
its  operational  characteristics. 

The  NOTS  Model  II  thrust  stand  was  designed  with  a  six-compoment 
capability  but  was  used  as  a  three -component  stand  for  this  development. 
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FIG.  2.  NOTS  Model  II  Thrust  Stand  With  Dummy  Load. 


An  evaluation  test  of  the  experimental  system  and  an  analog  simula¬ 
tion  study  have  been  conducted.  The  evaluation  test  is  discussed  on 
p age  35;  results  of  the  analog  simulation  study  are  given  in  Appendix 
A;  the  derivation  of  analog  simulation  equations  is  presented  In  Appen¬ 
dix  B;  Appendix  C  includes  the  derivation  of  system  transfer  functions. 


MASS  MEASUREMENT  TECHNIQUE 


The  mass  of  a  burning  rocket  motor  is  determined  by  continuously 
exciting  the  rocket-motor-stand  system  at  its  undamped  natural  frequency 
when  supported  on  springs  as  shown  schematically  in  Fig.  3*  If  the  vi¬ 
brating  system  is  assumed  to  have  a  single  degree  of  freedom  and  slowly 
varying  mass  then  the  mass  can  be  calculated  throughout  burning  from  the 
following  equation: 


M(t) 


K 

"n(t) 


(3) 


k 
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H4- 


FIG.  3* 


where 


=  undamped  natural  frequency  of  vibrating  system 


jj  /,  v  =  uuutuiijjeu  xiauux-fcu.  1 1  equexic.y  u 

'  (variable  with  burning  time) 


F(t)  =  FQsinci)t  (sinusoidal  force  generator),  lb 

M(t)  =  effective  vibrating  mass, 

K  =  lumped  spring  constant,  lb/ft 

lb-sec 


C  =  lumped  damping  constant, 


ft 


Z  =  vertical  displacement  of  rocket  motor  center  of 
gravity,  ft 


This  equation  is  derived  in  writing  the  dynamic  equation  of  motion. 
Applying  Newton fs  second  law  to  the  system  of  Fig.  3  and  assuming  a  non¬ 
changing  mass  environment,  the  differential  equation  of  translational 
motion  is  written: 

2 

F  sinwt  =  M  +  C  ^  +  KZ  (4) 

°  at2  dt 


Applying  the  LaPlace  transform  operator,  S 


-  sHs)  +  “(s)  *  E(s) 


+  w 


(5) 


5 


NAVWEPS  REPORT  S35U 


solving  for  Z^ 


J(s) 


F  a) 
o 


(S2  4-  0)2) 


(S  M  +  SC  +  K) 


rearranging 


J(S) 


F  w 

o 


M(S2  +  go2) 


<s+sH  +  S> 


(6) 


(7) 


By  definition,  any  linear  second-order  differential  equation  can  be 

written  in  terms  of  two  quantities:  t  and  go  . 

n 

The  damping  ratio  £  is  the  ratio  of  the  damping  that  exists  in  a 
second-order  system  to  the  critical  damping;  i.e.,  the  value  of  damping 
resulting  in  a  nonoscillatory  response  following  a  step  input.  The  un¬ 
damped  natural  frequency  u>n  is  the  frequency  of  oscillation  that  occurs 
for  zero  damping.  By  definition,  then 


2£u>  = 

n 


C 

M 


(8) 


=  § 
n  M 


(9) 


Rearranging,  Eq.  (9)  then  becomes  original  Eq.  (3): 
K 


M  = 


0) 


Substitution  of  Eqs.  (8)  and  (9)  in  Eq.  (7)  yields 


F  w 
o 


^S)  M(S2  +  <o2)  (S2  +  2£w  S  +  co2) 

n  n 


do) 


Rearranging 
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F  u) 
o 


J(S)  “  M(S2  +  a)2)  £(S  +  Swn)2  +  <4  (1  -  ?2)] 


(11) 


Using  the  appropriate  transform  function,  the  complete  solution  of 
the  displacement  of  the  rocket  motor  is  found  to  be: 


Z(t)  = 


F  ca 
o 


M(  cl+J-2w2w2-t4£2w2u)2 ) 2 
x  n  n  n  ' 


n 

-sin(cot-0  )  + 


-£u>  t 
n  .  / 
e  sin(co 


n 


2  (9  ' 

“V 


(12) 


where 


/  2£to  0)  \ 

(-a-S) 

\  CO  -  CO  y 


-2C  w2v/l  -l2 

2^2w2  -  w2  +  u>2 
n  n 


Considering  only  the  steady-state  solution  it  can  be  seen  that  as 
the  frequency  of  the  exciting  force  go  approaches  the  undamped  natural 
frequency  u>n,  the  phase  angle  approaches  90°  •  In  other  words,  when 
the  spring -mass  system  is  oscillating  at  the  undamped  natural  frequency 
the  displacement  Z  lags  the  exciting  force  in  phase  by  90 D  so  that 

for 


F(t)  -  FQsincot  (13) 

F 

Z(t)  =  - 2_  sin(cot  -  90°)  (l4) 

2M  £co 

n 


The  problem  remaining,  then,  is  to  design  a  control  system  to  main¬ 
tain  this  9°°  phase  relationship  continuously  throughout  burning.  By 
maintaining  this  relationship,  system  oscillation  at  the  undamped  natu¬ 
ral  frequency  is  insured  and  computations  of  mass  versus  burning  time 
can  be  made  by  using  Eq.  (3)* 


Gardner,  Murray  F.  and  John  L.  Barnes.  Transients  in  Linear  Sys¬ 
tems,  Vol.  1.  New  York,  Wiley,  1956.  Transform  pair  1.357;  Appendix  A. 
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EXPERIMENTAL  SYSTEM 


PHASE  DETERMINATION 

The  first  step  in  the  design  of  the  control  system  was  to  devise  a 
method  of  sensing  the  phase  relationship  of  the  sinusoidal -force  gener¬ 
ator  exciting  function  F(t)  with  that  of  the  spring- supported  rocket- 
motor  stand  system  displacement  function  Z(t). 

The  waveforms  in  Fig.  4a,  which  depict  the  phase  relationship  be¬ 
tween  these  two  functions  at  an  instant  in  time  when  the  spring -mass 


system  is  oscillating  just  below  the 
a  phase  angle  of  80°.  In  order 
to  detect  the  phase  polarity 
(i.e.,  lead  or  lag,  in  ref¬ 
erence  to  the  desired  9°°  value) 
a  third  waveform  is  added.  This 
waveform,  shown  in  Fig.  4b,  is 
always  leading  the  exciting  func¬ 
tion  by  90° .  For  the  condition 
shown  (i.e.,  for  the  system  os¬ 
cillating  below  the  undamped  nat¬ 
ural  frequency)  it  is  now  possible 
to  write  logic  equations  relating 
the  three  waveforms  and  the  10° 
phase  error. 

Using  Boolean  Algebra  termi¬ 
nology  and  selecting  A,__B,_C  to 
represent  positive  and  A,  B,  C  to 
represent  negative  values,  the 
following  equations  are  written: 


A  B  C  =  1 

(15) 

A  B  C  =  1 

(16) 

where  A  B  C  is  defined  as  A 
and  B  and  C. 

In  the  same  manner  it  is 
possible  to  write  the  logic 
equations  for  the  condition 
wherein  the  spring-mass  system 
is  oscillating  above  the  un¬ 
damped  natural  frequency.  As 
illustrated  in  Fig.  4c,  the 
logic  equations  are: 


undamped  natural  frequency,  show 


FIG.  4.  Phase  Comparator  Logic 
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A  B  C  =  1 

(17) 

A  B  C  =  1 

(18) 

These  equations  were  used  in  designing  the  diode  logic  circuitry 
(phase  comparator),  a  schematic  diagram  of  which  appears  in  Appendix  C* 
The  phase  comparator  generates  an  error  pulse  of  constant  amplitude  with 
the  width  equal  to  the  relative  phase  error.  The  pulse  is  positive  for 
a  system-oscillation  frequency  below  the  undamped  natural  frequency  and 
negative  for  frequencies  above  the  undamped  natural  frequency.  Since 
one  of  the  four  logic  equations  will  be  satisfied  at  each  crossover  of 
the  displacement  wave,  a  digital  sensor  exists  with  a  sampling  rate  of 
two  samples  per  cycle.  This  low  sampling  rate  in  a  closed-loop  system 
will  tend  to  increase  system-response  time  which  will  generally  degrade 
system  performance.  In  this  experimental  system,  however,  the  power 
train  has  a  time  constant  of  approximately  200  mill i sec  and  the  addi¬ 
tional  lag  introduced  by  the  lowest  sampling  rate  of  20  samples  per  sec¬ 
ond  does  not  appear  to  be  significant.  It  is  apparent  that  in  the  design 
of  fully  operational  systems  having  power  train-time  constants  consider¬ 
ably  lower  than  200  millisee,  techniques  for  increasing  the  digital  sam¬ 
pling  rate  will  have  to  be  devised. 


FORCE  GENERATOR 

The  force  generator  is  a  mechanical  device  consisting  of  two  sets 
of  eccentric  weights  connected  through  a  differential.  The  arrangement 
and  relative  rotation  of  the  weights  is  shown  in  Fig.  5* 


FIG.  5.  Rotating  Eccentric 
Weights.  Arrows  indicate 
direction  of  rotation. 


9 


NAVWEPS  REPORT  83$4 


With  an  individual  eccentric  weight  W^,  a  radius  of  rotation  rj_  and 
an  angular  velocity  w,  the  centrifugal  inertia  force  is  given  by 

W 

Fi  -  J  "S  (19) 

where 

F^  =  force,  lb 
Wp  =  effective  weight,  lb 
r-^  =  radius  of  rotation,  ft 
g  =  gravitational  force,  32.2  ft/sec2 
to  =  angular  velocity,  rad/sec 

Since  the  two  weights  of  each  set  are  coupled  together  and  rotate 
in  opposite  directions,  the  centrifugal  inertia  force  along  the  x  axis 
is  canceled  out  and  a  force  is  generated  only  along  the  y  axis.  The 
force  is  sinusoidal  and  has  a  maximum  magnitude  which  is  dependent  on 
the  phase  relationship  between  the  two  sets  of  weights.  Thus,  the  ex¬ 
citation  force  is  calculated  using  the  following  equation: 

4  w 

p  ss  - i  r  to2  sin(0n  -  eL)  sincot  (20) 

1  g  1  x  1  2'  x  J 

The  desired  magnitude  of  force  can  be  generated  by  selecting  the 
proper  phase  angle  between  the  two  sets  of  weights  with  the  differential; 
hence  the  term  sin^^-S^) . 

The  magnitude  of  the  force  necessary  to  excite  the  spring -mass  sys¬ 
tem  at  the  undamped  natural  frequency  is  dependent  upon:  (l)  the  accel¬ 
eration  limitations  in  the  vertical  plane  of  the  rocket  motor  under  test, 
(2)  the  damping  existing  in  the  spring-mass  system,  and  (3)  the  mass  of 
the  vibrating  system.  For  the  rocket  motor  and  stand  assembly  being 
tested,  the  following  specifications  apply: 
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Sprung  weight  with  live  rocket  motor 
(propellant  unburned)  in  stand 

Sprung  weight  with  rocket-motor  car¬ 
cass  (propellant  burned)  in  stand 

Effective  spring  constant 

Undamped  natural  frequency  (live 
motor  in  stand) 

Undamped  natural  frequency  (motor 
carcass  in  stand) 

Damping  ratio  £  (this  ratio  was  cal¬ 
culated  using  Eq.  (22)  and  esti¬ 
mating  F0,  while  maintaining  a 
1-g  acceleration) 

Maximum  vertical  acceleration 

Burning  time 

By  taking  the  second  derivative  of  Eq#  (l4),  the  following  expres¬ 
sion  is  obtained  for  the  acceleration  of  the  spring -mass  system  oscil¬ 
lating  at  the  undamped  natural  frequency: 

F 

*(t)  =  w cosut  (21) 

The  vibrating  structure  (including  a  first-stage  rocket  motor)  has 
a  mass  of  approximately  720  slugs  when  the  rocket  motor  propellant  is 
unburned  and  approximately  320  slugs  with  the  rocket  motor  carcass.  The 
damping  ratio  is  small  since  dash  pots  have  not  been  added  and  the  only 
damping  present  (such  as  propellant  flexibility)  is  that  which  is  inher¬ 
ent  in  the  system.  £  is  estimated  at  approximately  0.01.  Under  these 
conditions,  the  magnitude  of  the  exciting  force  can  be  computed  from 
Eq.  (21). 

Maximum  acceleration  occurs  when  coswt  =  1.  Therefore,  rewriting 
Eq.  (21) 


25,000  lb 

10,000  lb 
3-72X106  lb/ft 

10.88  cps* 

16.88  cps* 

0.01  or  1$ 

1.0  g 

60  sec 


Fq  =  2£MZ 


(22) 


Values  determined  experimentally. 
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and  for  the  conditions  above 

Fq  =  (2) (0. 01) (320) (32.2)  =  206  lb  (23) 

This  value  of  F0  will  result  in  a  1-g  acceleration  with  a  rocket 
motor  carcass  mounted  in  the  stand.  Since  £  is  considered  constant 
throughout  burning,  it  is  apparent  that  a  much  larger  exciting  force 
could  be  tolerated  under  live  motor  conditions.  However,  since  the 
force  generator  differential  will  not  be  operated  during  the  test,  it 
is  necessary  to  restrict  the  maximum  force  to  that  indicated  in  Eq.  (23). 

It  is  now  possible  to  calculate  the  maximum  horsepower  required  to 
drive  the  force  generator.  The  power  required  in  the  system  is  the  sum 
of  the  power  needed  to  provide  angular  acceleration;  to  overcome  damping 
in  the  spring -mass  system,  shaker,  and  drive  motor,  to  counteract  the 
losses  in  the  gear  train;  and  to  compensate  for  torque  feedback  due  to 
the  oscillatory  motion  of  the  shaker  in  relation  to  the  drive  motor. 

The  maximum  horsepower  requirement  is  calculated  by: 

T  =  JO  +  BQ  (unknown)  +  T«  +  gear  train  losses  (unknown) 

,  (2*0 

+  spring -mass  damping  v  ' 

and 

T f  =  mi  [£9  +  (Z  +  Z6T)  cos 0  +  ZB  sinS]  (25) 

where 

T  =  total  torque  reflected  into  the  drive  motor,  lb -ft 
T^  »  feedback  torque,  lb -ft 

0  =  angular  velocity,  rad/sec 

■  *  2 
G  =  angular  acceleration,  rad/sec 

2 

m£  =  effective  shaker  unbalance,  lb -sec 
Z  ~  effective  mass  lever  arm,  ft 

2 

J  =  effective  system  inertia,  lb -ft -sec 
The  derivation  of  Eq.  (25)  is  given  in  Appendix  B. 
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To  calculate  the  maximum  power  required,  9  and  9  will  be  assumed  to 
occur  simultaneously  and  have  the  following  values: 

9  max  =  106  rad/sec 

**  2 
9  max  =  10  rad/sec 


The  term  9  is  based  on  the  system  being  in  operation  with  a  rocket 
motor  carcass  in  the  stand  and  9  is  based  on  the  maximum  transient  re¬ 
sponse  requirement  as  indicated  in  the  figure  shown  on  page  30* 

Since  the  second  and  third  terms  of  Eq.  (25)  have  a  90°  relation¬ 
ship,  only  the  second  term  need  be  considered  to  obtain  peak  torque. 
When  cos#  =  1,  this  equation  becomes 

Tf  =  red  (i0  +  Z  +  ze2)  (2 6) 

where 

mi  =  0,02 
i  =  0.096 

then 


Tf  =  (0.02) (O.O96  x  10  +  32.2  +  0.00287  x  1062)  =  2.60  lb -ft 
The  additional  torque  required  to  provide  the  desired  acceleration  is 
J 9  =  (10)(0.71)  =  7-1  lb -ft 


where 


J  =  0.71 


Then 


T  =  +  J9 =  2.60  +  7*1  =  9*7  lb-ft  total  known  torque 
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The  horsepower  transmitted  by  a  shaft  making  n  revolutions  per 
minute  under  a  torque  of  T  lb-ft  is 


2nn  T 

~  33,000 

(27) 

=  (2>t)(l6.88)(6o)(9.7)  _  i  ft7 

q "3  nnn  x  *u  > 


The  horsepower  requirement  to  provide  acceleration  and  to  compensate  for 
torque  feedback  then,  is  I.87.  The  additional  power  requirement  to  over¬ 
come  spring -mass  damping  is  obtained  from  the  following  power  equation, 
where  Z  is  the  velocity  of  the  displacement  Z  and  C  is  the  damping  con¬ 
stant  as  in  Eq.  (2): 


Power  =  Z  C 


(28) 


The  damping  C  as  determined  from  Eq.  (8)  is  680 
placement  velocity  Z  is  obtained  by  differentiating  Eq,  (l4);  therefore, 

Power  =  (0.30^)2(680)  =  63  — —  — 


HP  = 


il 

550 


0.12 


This  results  in  a  total  known  instantaneous  peak  power  requirement 
of  approximately  two  horsepower.  A  drive  motor  is  normally  selected  by 
considering  both  the  instantaneous  and  average  power  requirements.  In 
this  case,  however,  the  average  power  dissipated,  other  than  that  due  to 
damping,  is  difficult  to  ascertain  analytically.  It  has  been  determined 
experimentally  that  the  average  power  dissipated  is  approximately  one 
horsepower.  It  appears,  therefore,  that  the  selection  of  the  5-horse¬ 
power  motor  used  in  this  experiment  is  somewhat  conservative  and  that 
improved  response  could  be  realized  by  carefully  selecting  a  motor  just 
matching  the  power  requirements. 
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SYSTEM  ANALYSIS 


The  block  diagram  of  Fig.  6  represents  the  complete  experimental 
control  system  when  it  is  operating  as  a  relay  type  with  sampling*  The 
system  was  initially  constructed  as  a  straight  sampling  type  (i.e.,  as 
shown  in  Fig.  6,  without  the  flip-flop).  Using  the  sampled  data  system, 
a  rocket-motor  firing  was  conducted;  the  results  are  discussed  later  in 
this  report.  The  system  performed  satisfactorily;  however,  in  prepara¬ 
tion  for  a  second  evaluation  test  it  was  found  that  the  addition  of  the 
flip-flop  improved  system  perf ormance .  (it  was  planned  to  conduct  the 
second  test  using  the  configuration  shown  in  Fig.  6;  however,  this  test 
could  not  be  conducted  because  of  a  rocket  motor  malfunction  which  par¬ 
tially  destroyed  the  drive  motor  and  shaker.) 


SINE-COSINE 

POTENTIOMETER 


FIG.  6.  Block  Diagram  of  the  Experimentally  Tuned  Control  System. 


Although  the  evaluation  test  was  conducted  with  the  system  operat¬ 
ing  as  a  straight  sampling  type,  it  is  felt  that  significant  improvement 
in  closed-loop  performance  was  obtained  with  the  addition  of  the  flip- 
flop  and  that  the  analysis  should  be  made  for  the  relay  sampling  type 
system.  No  recorded  data  of  the  relay- sampling  system  were  obtained  be¬ 
cause  of  the  rocket-motor  malfunction.  The  system  was  operated,  however, 
with  a  rocket -motor  carcass  in  the  stand  prior  to  this  catastrophe  and 
the  system  performance  was  observed.  The  significant  performance  char¬ 
acteristics  are  indicated  for  comparison  on  the  response  curves  obtained 
with  the  straight  sampling  system. 


15 


NAVWEPS  REPORT  835*+ 


In  Fig.  6 ,  comparator  No.  1  and  the  ’sampler1  represent  the  phase 
comparator  described  previously.  The  phase  of  the  waveform  0,  represent¬ 
ing  the  exciting  force  is  compared  with  the  phase  of  the  waveform  Z  rep¬ 
resenting  the  displacement  of  the  oscillating  spring-mass  system.  A 
sampling  rate  of  two  samples  per  cycle  has  been  established.  The  output 
of  the  sampler  is  used  to  trigger  a  flip-flop  which  in  turn  is  used  as 
the  error  signal  for  integration.  (At  least  one  stage  of  integration 
in  addition  to  that  performed  by  the  drive  motor  is  desirable  in  veloc¬ 
ity  control  since  this  results  in  a  Type  2  system  which  exhibits  a  zero 
steady-state  error  in  a  simple  configuration — i.e.,  a  single  loop  with 
continuous  feedback.  A  theoretical  derivation  is  given  in  Chapter  7  of 
Servomechanism  Analysis.*)  The  addition  of  the  flip-flop  changes  the 
servomechanism  from  a  straight  sampling  type  to  a  relay  type  with  sam¬ 
pling. 

Usually,  the  performance  characteristics  of  relay  servos  are  not 
particularly  favorable;  that  is,  a  servo  with  continuous  (proportional) 
feedback  will  generally  out-perform  the  relay  type.  In  this  case,  how¬ 
ever,  where  it  is  necessary  to  compare  the  phase  relationship  of  two 
sinusoidal  waveforms  (one  of  which  has  a  variable  amplitude),  a  practi¬ 
cal  device  for  comparing  the  phase  relationship  continuously  is  not 
readily  available.  The  phase  comparator  thus  uses  diode  logic  circuitry 
to  produce  pulses  which  are  phase  dependent  at  each  zero  crossover  of 
the  displacement  waveform.  The  pulses  are  of  constant  amplitude  but 
have  a  width  directly  proportional  to  phase  error. 

Comparing  the  system  performance  during  experimental  operation, 
both  as  a  straight  sampling  type  and  as  a  relay  type  with  sampling,  it 
was  found  that  the  relay  type  exhibited  superior  characteristics.  With 
the  relay  configuration,  it  was  possible  to  more-than-double  the  forward- 
loop  gain  of  the  control  system  with  a  consequent  improvement  in  fre¬ 
quency  response.  The  steady-state  oscillation  or  limit  cycle  of  the 
servo  system  was  more  pronounced  when  operated  as  a  relay  type;  however, 
the  frequency  of  oscillation  was  higher,  and  better  response  was  obtained 
for  small  step  inputs  in  phase  error.  It  should  be  noted,  however,  that 
with  the  use  of  highly  sophisticated  compensating  circuitry,  a  straight 
sampling  system  performance,  comparable  to  that  of  the  relay  type,  could 
very  likely  be  obtained.  The  compensating  circuitry  necessary  to  attain 
this  result  is  quite  complex  making  experimental  improvement  of  closed- 
loop  system  performance  extremely  difficult.  While  on  the  other  hand, 
the  relay  system  required  relatively  elementary  compensating  circuitry 
and  produced  the  desired  result. 

From  the  control  system  diagrammed  in  Fig.  6,  it  is  apparent  that 
several  nonlinearities  exist.  (A  component  is  said  to  be  nonlinear  if 
the  equations  describing  its  dynamics  contain  terms  made  up  of  the  pro¬ 
duct  or  quotient  of  a  dependent  variable  and  a  time  variant  coefficient.) 
The  equation  describing  the  spring-mass  dynamics  is  written: 

'  ‘ 

Thaler,  G.  J.  and  R.  G.  Brown.  Servomechanism  Analysis.  New  York, 
McGraw-Hill,  1953* 
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F  (^)2 
0  W 


sinU(t)t  = 


M 


(t) 


a2z 

at2 


P  dZ 

0  at 


+  KZ 


(29) 


Since  this  equation  is  of  the  nonlinear  type,  it  cannot  be  handled 
by  regular  LaPlace  methods. 

The  complete  solution  of  the  nonlinear  system  will  not  be  attempted, 
only  the  two  end  points  of  operation  will  be  considered.  That  is,  the 
control  system  will  be  analyzed  operating  first  in  a  nonchanging  mass 
environment  with  a  rocket  motor  carcass  and  then  with  a  live  rocket  mo¬ 
tor  mounted  in  the  thrust  stand.  To  simplify  and  linearize  the  system, 
the  spring-mass  dynamics  will  be  considered  analogous  to  a  low-pass  fil¬ 
ter.  That  is,  the  spring -mass  system  will  be  considered  to  introduce 
phase  lag  and  amplitude  attenuation  in  the  feedback  loop.  Data  obtained 
with  the  experimental  system  and  in  the  analog  simulation  studies  (Ap¬ 
pendix  A),  indicate  that  a  lag  in  the  spring -mass  system  to  simulated 
step  inputs  exists  and  that  the  response  curve  approaches  that  of  an 
exponential . 

This  lag  characteristic  is  described  in  some  detail  in  NAVWEPS 
Report  7741*  which  illustrates  that  a  step  change  in  mass  will  result  in 
an  exponentially  rising  input  error  to  the  control  system.  The  assump¬ 
tion  has  been  made  here  that  this  lag  will  also  appear  when  changes  in 
driving -function  frequency  are  experienced  and  that  this  lag  is  in  the 
feedback  path.  This  assumption  is  open  to  considerable  discussion;  how¬ 
ever,  it  is  shown  that,  with  this  assumption,  reasonable  agreement  with 
experimental  data  can  be  obtained.  Using  the  experimentally  determined 
time  constants,  the  following  transfer  functions  will  be  used  to  simu¬ 
late  the  spring-mass  system  in  the  analysis  of  the  two  end  points  of 
operation: 


with  a  rocket  motor  carcass  in  the  stand 

z  -  3 

S  S  +  3 

with  a  live  rocket  motor  in  the  stand 

Z  _  2 

e  s  +  2 

This  simplifies  the  system  block  diagram  to  that  of  Fig.  J. 

•x-  - 

U.  S.  Naval  Ordnance  Test  Station.  Determination  of  Rocket-Motor 
Mass  by  Measurement  of  the  Natural  Frequency  of  a  Mass-Spring  System,  by 
Benjamin  Glatt.  China  Lake,  Calif.,  NOTS,  15  June  1961.  (NAVWEPS  Re¬ 
port  7741,  NOTS  TP  2706) . 


(30) 


(31) 
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FLIP-  DRIVE 

SAMPLER  FLOP  INTEGRATOR  COMPENSATION  PRE-AMP  MAG-AMP  MOTOR 


FIG.  7-  Simplified  Block  Diagram  of  the  Compensated 
Control  System. 


It  should  be  pointed  out  that  these  time  constants  may  be  in  con¬ 
siderable  error.  As  noted  in  NAVWEPS  Report  the  time  constants 

are  dependent  on  the  damping  and  mass  of  the  system  at  the  time  the  step 
change  in  mass  is  initiated.  The  experimental  data  were  obtained  by 
suddenly  changing  the  driving-function  frequency  and  observing  the  lag 
of  the  spring-mass  system.  Since,  in  this  experiment,  both  the  lag  of 
the  power  train  and  the  spring -mass  system  are  combined,  some  error  is 
evident  in  the  interpolation.  It  is  also  noted  that,  although  the  time 
constant  determined  under  live-motor  conditions  compares  favorably  with 
that  of  the  analog  study,  there  is  a  considerable  difference  between  the 
time  constant  obtained  under  rocket-motor  carcass  conditions  and  that 
obtained  in  the  analog  study. 

In  the  analog  simulation,  it  was  necessary  to  vary  the  damping  with 
burning  time  and  to  distribute  it  unequally  between  front  and  rear  sup¬ 
ports  to  provide  the  needed  coupling  to  cause  pitching  and  thus  repro¬ 
duce  the  displacement  curves  obtained  experimentally.  It  has  since  been 
determined  under  system  operation  that  the  structural  oscillation  of  the 
thrust  take-out  member  was  also  introducing  pitching.  It  is  evident, 
therefore,  that  the  damping  used  in  the  analog  study  for  the  rocket - 
motor  carcass  condition  is  in  error  and  that  the  time  constant  of  the 
rocket-motor  lag  for  this  condition  is  also  in  error. 

The  speed  of  the  shaker -drive  motor  is  controlled  by  regulating  the 
armature  current  while  maintaining  a  constant  field  strength.  The  dy¬ 
namic  equations  pertaining  to  this  type  of  control  are  well  developed.* 
The  transfer  function  of  the  motor  when  driving  an  inertia  load  and  neg¬ 
lecting  frictional  losses  is: 


* 

Thaler,  G.  J.  and  R.  G.  Brown. 


Op.  cit.,  pp.  62-65. 
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K 

JR 


S  S  + 


jr 

JR 


(32) 


where 


K  =  motor  constant, 


ohm-sec 

rad 


R  =  motor  armature  resistance,  ohms 
a  9 

2 

J  =  total  inertia  load,  ft -lb -sec 

For  the  system  under  development,  the  motor  constants  K,  and  Ra 
were  determined  experimentally  (Appendix  C).  The  transfer  function  is: 


G 


e6  S(S  +  10) 


(33) 


The  magnetic  amplifier  dynamic  characteristics  were  also  determined 
experimentally  (Appendix  C)  and  the  transfer  function  is: 


e6  A2 

i5  (S  +  10) 


(34) 


Tachometer  feedback  is  utilized  in  a  minor  loop*  Since  this  is  a 
velocity  servo,  the  derivative  of  the  tach  or  9  is  actually  fed  back 
and,  since  the  tach  output  contains  high-frequency  components  in  the 
form  of  ripple,  a  filter  is  employed  with  the  differentiator.  The  fil¬ 
ter  is  connected  in  the  feedback  circuit  of  the  operational  amplifier 
used  as  the  differentiator;  filtering  is  accomplished  since  the  gain  be¬ 
comes  a  function  of  frequency.  The  lead  circuit  is  used  to  reshape  the 
filter  characteristics  slightly  and  to  allow  optimization  of  the  tach 
circuit  under  operating  conditions.  Derivation  of  the  transfer  functions 
appear  in  Appendix  C . 

The  steady -state  transfer  function  of  the  mag -amp  and  drive  motor 
as  a  unit  was  obtained  experimentally  to  determine  if  there  were  any 
nonlinearities  in  this  portion  of  the  control  system.  The  curve  plotted 
in  Fig.  8  shows  the  nonlinearity  that  is  present.  If  the  mag -amp  and 
drive-motor  transient  characteristics  are  considered  not  to  change  from 
the  experimentally  obtained  values  then  the  analysis  of  the  end  points 
of  operation  can  continue  by  assuming  an  appropriate  change  in  the  for¬ 
ward  loop  gain  as  dictated  by  the  curve  in  Fig.  8. 
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Q  (CPS) 


FIG.  8.  Mag -Amp  and  Drive  Motor  Static  Characteristics. 


Under  closed -loop  operation  it  was  oh served  that  the  control  system 
exhibited  a  limit-cycle  frequency  of  approximately  one  cycle  per  second. 
The  sampling  frequency  at  the  low  frequency  end  (with  a  live  rocket  mo¬ 
tor  mounted  in  the  stand)  is  about  20  samples  per  second  and  about  36 
samples  per  second  at  the  high  frequency  end.  Since  the  sampling  rate 
is  at  least  twenty  times  the  switching  rate  of  the  relay  (flip-flop),  in 
this  analysis  the  sampler  will  be  considered  to  have  an  infinite  fre¬ 
quency*  That  is,  the  error  detector  will  be  considered  to  be  continuous. 

The  describing  function  technique  for  analysis  of  relay  servos  is 
particularly  useful  when  relays  approaching  the  ideal  are  used  and  when 
the  system  is  of  a  high  order.*  A  linear  approximation  is  made  in  the 
frequency  domain  by  defining  a  describing  function  in  terms  of  the  Fourier 
series  for  the  component  response  to  a  sinusoidal  input.  With  an  ideal 
relay,  the  application  of  a  sinusoidal  input  will  result,  in  a  Fourier 
series  describing  the  output  consisting  of  the  fundamental  input  and  all 
odd  harmonics  thereof.  The  approximation  is  made  with  the  assumption  that 
the  relay  output  consists  only  of  the  fundamental  sinusoidal  input.  For 
this  to  hold,  it  is  necessary  that  the  feedback  signal  generated  by  the 
relay  output  be  a  pure  sinusoid.  Since  a  relatively  high-order  system 
exists  (i.e.,  one  with  considerable  low-pass  filtering),  this  assumption 
appears  valid.  Since  the  ideal  relay  presents  zero  phase  shift,  the  fol¬ 
lowing  described  transfer  function  is  defined  in  amplitude  only: 


* 

Thaler,  G.  J.  and  Marvin  P.  Pastel.  Analysis  and  Design  of  Non¬ 
linear  Feedback  Control  Systems.  Chapter  4.  New  York,  McGraw-Hill,  1962. 
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Describing  function 


Vu) 

A 


L* 


where 


(35) 


A  =  amplitude  of  input  sinusoid  to  cause  saturation 


F(A,u>) 


=  amplitude  of  output  sinusoid 


0  as  phase  angle  "between  input  and  output  functions 


For  the  system  under  study, 
of  0.1  volt  and  an  output  F^ 


assume  a  maximum  input  sensitivity  A 
of  30  volts.  Then 


=  300  [_0^ 


(36) 


The  block  diagram  of  the  servo  system  is  now  reduced  to  that  shown 
in  Fig.  9*  It  is  observed  that  the  effect  of  the  nonlinear  component 
Gj)  on  the  control  system  is  to  introduce  a  variable  forward  loop  gain. 
When  large  errors  and  transients  exist  in  the  control  system  the  forward 
loop  gain  is  minimized  since  the  input  A  of  the  describing  function, 

Eq.  (35),  increases  while  the  output  F  remains  constant.  The  maximum 
forward  loop  gain  occurs  under  steady-state  operation  when  the  control- 
system  feedback  is  just  sufficient  to  trigger  the  relay  (flip-flop) . 


DESCRIBING  INTEGRATOR  a 

FUNCTION  COMPENSATION 


LIVE  ROCKET  MOTOR 
IN  STAND 


ROCKET  MOTOR 
CARCASS  IN  STAND 


\/f-  2 


i/r  *  3 


A?K2  =  1190 


A2K2=2580 


FIG.  9*  Diagram  Showing  Further  Simplification  of  the 
Compensated  Control  System. 
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ROOT  LOCUS  STABILITY  ANALYSIS 

Once  the  describing  function  has  been  defined  it  is  possible  to 
examine  the  servo  performance  by  the  root  locus  method.  Since  the  non¬ 
linear  element  is  considered  to  introduce  variable  gain  only,  the  locus 
of  roots  will  not  be  affected  and  only  a  single  overall  curve  is  required 
to  depict  servo  performance  as  a  function  of  forward  loop  gain. 

The  uncompensated  control  system  of  Fig*  10  yields  the  root-locus 
plots  of  Figs.  11  and  12.  It  is  evident  that  the  uncompensated  system 
is  unstable  for  all  values  of  gain.  It  is  also  evident  that  the  addi¬ 
tion  of  zeros  along  the  negative  real  axis  (i.e.,  (S  +  A)  terms  in  the 
numerator  of  KGHp)  will  tend  to  shift  the  root  locus  into  the  left-half 
plane  and  in  general  stabilize  the  system.  Several  compensating  circuits 
are  available  which  can  be  used  in  the  forward  loop  to  obtain  these  zeros. 
The  final  compensation  of  the  control  system,  as  indicated  in  Fig.  9}  was 
determined  experimentally.  It  is  noted  that  the  addition  of  the  tach 
feedback  circuit  contributes  two  additional  zeros  in  the  open  loop  trans¬ 
fer  function.  This  is  indicated  in  the  equivalent  series  element  in 
Fig.  13. 


FIG.  10.  Block  Diagram  of  the  Uncompensated 
Control  System. 


In  the  analysis  of  the  compensated  system  of  Fig.  9  the  root-locus 
curve  of  the  minor  loop  is  plotted  first  and  the  roots  of  the  system  for 
a  particular  minor  loop  gain  are  determined.  In  turn,  these  roots  are 
used  as  poles  to  write  the  overall  open -loop  transfer  function.  The  root 
locus  of  the  complete  system  is  then  plotted  and  the  point  of  system  oper¬ 
ation  is  determined. 
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FIG.  12.  Root  Locus  of  the  Uncompensated  System  With 
a  Live  Rocket  Motor  in  the  Stand. 
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DESCRIBING 

FUNCTION 


i/r  =  2 

A2  K2  =  ||90 
S,  s  —10.2  +  j  4.0 
S- =  — 10.2  —  j  4.0 
Sj=-I.7  4-  j  0.6 
S4  =  —1.7  —  j  0.6 


\/r*  3 

A2K2=25B0 
S|  =  — 10,4  +j  6.2 
S2  =  -I0.4~j  6.2 
S3=  “  1.6  -+•  J  0.7 
S4=~  1.6-  J0.7 


FIG.  13*  Final  Block  Diagram  of  the  Compensated  Control  System 
With  the  Equivalent  Series  Transfer  Function  Replacing  the  Minor 
Loop. 


The  system  will  be  analyzed  first  using  the  system  constants  asso¬ 
ciated  with  a  rocket  motor  carcass.  The  minor  open-loop  transfer  func¬ 
tion  of  Fig.  9  is  given  as: 


kgh2 


28  S(S  +  1) 

(S  +  10)2(S  +  2)2 


(37) 


A2K2K4  =  28 


The  roots  at  the  point  of  operation  as  determined  from  the  locus  of 
roots  shown  in  Fig.  l4  are  as  follows: 


=  -10.4  +  j6. 2 
S2  =  -10.4  -  j6. 2 
~  "  1.6  +  jO.Y 
S4  =  -  1.6  -  j0.7 


The  minor -loop  transfer  function  can  now  be  written  as  follows: 


Q  A2K2  (S  +  2)2 

\  =  S(S  +  S^S  +  S2)(S  +  S3)(S  +  Sk) 


(38) 
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FIG.  l4.  Root  Locus  of  Minor  Loop. 


This  transfer  function  can  now  he  handled  as  a  series  element  in 
the  outer  loop  as  illustrated  in  Fig.  13*  The  outer  open-loop  transfer 
function  is  written: 


0.9  A  K  (S+l)2(S+2)2 

KGH  =  - 5 - <=-* - - 

7  Sc:(S+ll)(S+50)(S4i)(S451)(S+S2)(S+S3)(S+S4) 

For  the  motor  carcass  condition  =  2,5^0  and  7  -  1/3  seconds. 

KGH.  -  - KSP  !**&*£)! -  (HO) 

s  ( s+50 ) ( s+3 ) ( s+11 ) ( s+s1 ) ( S+S2 ) ( s+s3 ) ( S+S4 ) 

The  root-locus  plot  for  this  open-loop  transfer  i  unoKon  is  shown  in 
Fig.  15. 

The  same  procedure  is  followed  in  plotting  the  overall  root  locus 
with  a  live  rocket  motor  in  the  stand.  Since  the  minor -loop  gain  is 
changed  for  this  condition,  a  new  KGH2  must  he  plotted  and  another  set 
of  roots  determined.  Thus, 
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FIG.  15.  Root  Locus  of  the  Compensated  System  With  a 
Rocket  Motor  Carcass  in  the  Stand. 


kgh2 


12.9  s(s  +  1) 
(S  +  2)2(s  +  10)2 


(41) 


a2k2k4  =  12.9 

The  roots  at  the  point  of  operation  as  determined  from  the  locus  of  roots 
in  Fig.  l4  are  as  follows: 


S  =  -10.2  +  j4.0 

52  =  -10.2  -  j4.0 

53  «  -  1.7  +  JO. 6 
Sk  «  -  1.7  -  JO. 6 


The  minor -loop  transfer  function  can  now  he  written: 


This  transfer  function  is  now  handled  as  a  series  element  in  the 
outer  loop  and  the  following  open-loop  transfer  function  is  written: 
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KGH1  = 


3340  (S-t-l)2  (S+2)2 


S  (S+ll)(S+50)(S+2) (S+S^) (S+S2) (S+S3) (S+S^) 


(43) 


The  root-locus  plot  for  this  open-loop  transfer  function  is  shown 
in  Fig,  1 6.  It  is  interesting  to  note  that  the  loci  of  Figs.  15  and  1 6 
are  dependent  only  on  the  linear  system.  The  effect  of  the  variable - 
gain  nonlinearity  (describing  function)  is  to  move  the  root  position 
along  the  locus. 


FIG.  1 6.  Root  Locus  of  the  Compensated  System 
With  a  Live  Rocket  Motor  in  the  Stand. 


In  a  relay  servo  system  of  the  type  under  analysis,  maximum  gain 
occurs  under  minimum  error  conditions.  It  is  seen  that  the  roots  for 
this  condition  are  located  at  their  greatest  displacement  from  the  open- 
loop  poles.  For  a  large  Initial  disturbance,  such  as  a  step  change  in 
mass,  the  gain  is  low  and  the  system  will  exhibit  a  low  natural  frequency. 
It  is  evident,  then,  that  in  this  type  of  servo  system  the  maximum  gain 
operating  point  should  produce  a  system  damping  ratio  somewhat  lower  than 
the  nominal  value  of  0.7*  The  system  under  study  was  experimentally 
tuned  at  the  high-frequency  end  point  (or  in  accordance  with  the  root 
locus  of  Fig.  15)*  The  approximate  operating  point  indicates  a  damping 
ratio  of  0.44  and  an  undamped  natural  frequency  of  5*5  rad/sec. 
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It  should  he  pointed  out  that  considerable  error  may  be  present  in 
the  calculation  of  the  maximum  describing  function  gain.  It  is  felt 
that  the  experimentally  determined  threshold  sensitivity  of  the  phase- 
comparator  flip-flop  configuration  may  be  in  error  by  as  much  as  50$  • 

An  increase  in  forward-loop  gain  of  this  magnitude  would  greatly  in¬ 
crease  the  natural  frequency  but  would  tend  to  make  the  system  unstable. 
An  interesting  point  in  favor  of  this  type  of  servo  system,  however,  is 
its  ability  to  operate  very  near  the  imaginary  axis  while  maintaining 
complete  stability.  This  is  best  understood  by  visualizing  the  large 
change  in  gain  which  takes  place  when  the  system  is  oscillating  about r 
the  zero  error  point.  At  the  minimum  error  pbint  the  describing  func¬ 
tion  gain  is  approximately  300.  This  corresponds  to  a  phase -angle  error 
(i.e.,  phase  angle  difference  from  90°  between  the  forcing  function  and 
displacement  function)  of  about  1°.  Under  maximum  error  conditions,  or 
at  a  1$  frequency  error,  a  phase-angle  error  of  about  10°  is  evident. 
This  reduces  the  describing  function  gain  to  30.  Thus,  a  total  forward- 
loop  gain  reduction  by  a  factor  of  10  has  taken  place.  This  tends  to 
stabilize  the  system  by  sluggishly  reducing  the  phase  error  to  zero. 


PERFORMANCE  CHARACTERISTICS 

Because  of  the  rocket-motor  malfunction  that  precluded  completion 
of  the  second  evaluation  test,  limited  data  were  obtained  with  the  sys¬ 
tem  operating  as  a  relay  type  with  sampling.  The  data  obtained  were 
observed,  but  not  recorded.  The  steady-state  performance  characteris¬ 
tics  were  observed  using  a  rocket-motor  carcass  as  the  mass  environment. 
The  experimental  results  are  interpolated  on  Fig.  17*  Figures  17  and  18 
also  show  curves  of  the  experimental  results  obtained  with  the  system 
operating  as  used  in  the  first  evaluation  test;  i.e.,  as  a  straight  sam¬ 
pling  type. 

To  obtain  the  curves  of  these  figures,  the  control  system  (straight 
sampling  type)  was  first  operated  using  a  live  rocket  motor  as  the  mass 
environment.  A  step  change  in  mass  was  approximated  by  closing  the  for¬ 
ward  loop  as  the  system  was  oscillating  out  of  resonance  with  phase  er¬ 
ror.  The  system  responded  to  peak  overshoot  in  approximately  1.4  sec 
and  settled  to  within  1$  of  the  undamped  natural  frequency  in  about  3 
sec.  The  driving -function  response  for  this  condition  is  plotted  in 
Fig.  18.  A  steady-state  oscillation  about  the  natural  frequency  of  l/3 
cps  was  observed  with  an  0.8$  maximum  deviation  from  the  mean. 

With  a  rocket-motor  carcass  mounted  in  the  stand,  the  system  re¬ 
sponded  to  peak  overshoot  in  approximately  1*5  sec  and  settled  to  within 
1$  of  the  undamped  natural  frequency  in  approximately  2.1  sec.  A  steady- 
state  oscillation  about  the  natural  frequency  of  l/2  cps  was  observed 
with  a  0.5 $  maximum  deviation  from  the  mean.  The  driving -function  re¬ 
sponse  for  this  condition  is  shown  in  Fig.  17* 
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It  is  evident  from  the  curves  of  Fig,  17  that  the  system  exhibits 
improved  transient : characteristics  when  operated  as  a  relay  type  with 
sampling  and  maintains  steady-state  characteristics  comparable  to  those 
observed  when  the  system  was  operated  as  a  straight  sampling  type.  The 
root  locus  plot  of  Fig.  1 6  indicates  that  the  transient-response  charac¬ 
teristics  will  be  somewhat  less  desirable  with  a  live  rocket  motor  in 
the  stand  than  it  was  with  the  rocket -motor  carcass  due  to  the  lower 
gain  and  resonant  frequency  of  the  system  at  the  operating  point. 

Techniques  have  been  developed  for  theoretically  approximating  the 
limit-cycle  frequency  and  transient  response  of  relay  servos,*  The  tech¬ 
nique  described  in  the  Analysis  and  Design  of  Nonlinear  Feedback  Control 
Systems  incorporates  the  use  of  the  describing  function  and  root -locus 
plot.  An  iterative  process  is  involved  and  although  some  degree  of  ac¬ 
curacy  is  attainable,  the  analysis  becomes  very  laborious  in  high-order 
systems.  In  general,  for  practical  purposes,  the  techniques  which  have 
been  developed  for  transient  analysis  of  nonlinear  servos  are  limited  to 
third-order  systems  or  lower.  As  a  very  rough  approximation,  however, 
it  is  possible  to  make  the  following  calculation.  Consider  the  equation 
given  for  the  resonant  frequency  of  a  lightly  damped  second -order  system: 


wr  =  wn'\/1  "  ^  (44) 

In  a  lightly  damped  system  the  transient  time  from  t  =  0  to  peak  over¬ 
shoot  is  approximately  equal  to  one -half  cycle  of  the  resonant  frequency 
oscillation  time.  Therefore,  from  the  root  locus  plot  in  Fig.  15: 


wr  =  5.5J1-0.44" 


(45) 


From  this 


t  =0.7  sec 

P  1 


From  Fig.  1 6: 


wr  =  3.6J1-0.55  ,  and 


(46) 


t  =  1.25  sec 
P 


Thaler,  G.  J.  and  Marvin  P.  Pastel.  Op.  cit.,  pp.  183-190. 
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FIG.  17*  Experimental  System  Response 
as  a  Function  of  Time  to  a  Simulated 
Step  Input  (Rocket  Motor  Carcass) . 
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FIG.  18.  Experimental  System  Response  as  a 
Function  of  Time  to  a  Simulated  Step  Input 
(Live  Rocket  Motor) ,  System  Operating  as  a 
Straight  Sampling  Type. 
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A  measure  of  agreement  is  evident  between  the  value  of  0*7  seconds  and 
the  experimental  data.  It  is  reasonable  to  assume  that  the  value  of 
1.25  seconds  is  equally  representable  and  that  the  transient  response 
of  the  relay  type  system  would  also  show  improvement  over  the  straight 
sampling  type  at  the  low  frequency  operating  point* 

In  developing  a  high-order  system  with  relay  and  other  nonlinear¬ 
ities,  theoretical  calculations  must  be  used  only  as  a  building  block. 
That  is,  the  theory  is  used  to  establish  a  first  approximation  from  which 
components  having  a  wide  controllable  range  of  values  are  designed  or  se¬ 
lected.  In  system  optimization,  the  root-locus  plots  become  very  valu¬ 
able  in  pointing  out  the  compensation  needed  for  improved  performance. 

For  instance,  several  possibilities  for  further  system  refinement  can  be 
listed  from  the  root  locus  plots  of  Figs.  15  and  l6.  The  difference  be¬ 
tween  these  plots  can  be  attributed  mainly  to  the  change  in  minor -loop 
gain  in  traversing  from  one  operating  point  to  the  other.  It  follows, 
therefore,  that  an  increase  in  the  minor -loop  gain  would  very  probably 
improve  system  performance  even  further.  The  primary  restriction  to  this 
is  steady-state  error  and  instability  both  of  which  are  affected  by  for¬ 
ward-loop  gain.  The  main  objective  is  to  establish  the  maximum  undamped 
natural  frequency  for  an  allowable  steady-state  error  thus  optimizing  the 
system rs  transient  response.  Modification  of  the  compensating  circuit  in 
the  forward  loop  could  be  made  to  shift  the  root  locus  more  into  the  neg¬ 
ative  plane.  This  would  also  Improve  the  transient  response. 

Additional  experimental  data  were  obtained  to  determine  the  validity 
of  the  assumption  that  a  single  degree  of  freedom  system  existed. 

Figure  19  shows  a  plot  of  the  front  and  rear  displacements  of  the 
stand  legs  into  the  steel  springs  with  a  rocket  motor  carcass  in  the 
stand;  Fig.  20  plots  the  displacement  with  a  live  rocket  motor  in  the 
stand.  Although  the  stand  was  statically  decoupled,  it  is  apparent  that 
a  two-degree-of -freedom  system  existed  consisting  of  two  modes — pitching 
and  translational.  A  look  at  the  equations  for  this  type  of  system  is 
revealing.  The  equations  for  a  completely  coupled  two -degree -of -freedom 
system  (Appendix  B)  are  as  follows: 

M(t)Z  +  (Cf4Cr)  Z  +(£rCr-£fCf)f  +  (Kf+K.)Z  +  (yg-ifKf)t 

p 

=  sino)  t  (47) 

J(t)f  +  (&«&)*  +  (-ercr-*fcf)z  +  (^Kf+/Kr)f 

+  (JrWr )z  =  0  (48) 
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By  statically  decoupling  the  system  (i.e.,  by  adjusting  K  and  Kf  so 

that  K  =0),  the  equations  are  simplified  to  Eqs.  (4-9)  and  (50). 

1  x  r  r 

The  center-of -gravity  shift  was  found  experimentally  to  be  small  so  that 

Z  and  Z-  remain  constant.  Therefore, 
r  f  7 

M(t)Z+(Cf+Cr)Z+(-erCr-ifCf)t  +  (Kf+Kr)Z  =  mA^sin^t  (49) 

j(t  )tMijCf+J®cr)t  +  (ircr-ifcf)z  +  (i^cf+Ar)t  =  0  (50) 


This  statically  decoupled  system  is  still  coupled  dynamically 
through  damping,  resulting  in  the  data  shown  in  Figs.  19  and  20.  It 
was  later  determined  that  the  thrust  takeout  member  was  also  contribut¬ 
ing  a  pitching  moment. 

The  major  trouble  caused  by  the  pitching  mode  is  that  at  some  time 
during  burning  the  pitching-mode  frequency  is  equal  to  the  translational 
mode.  At  this  point,  phase  errors  are  introduced  into  the  translational 
sensors  and  the  control  loop  will  not  function  properly.  Since  it  was 
not  possible  during  design  to  compute  the  actual  damping  encountered, 
Figs.  19  and  20  represented  the  first  concrete  information  obtained  on 
damping.  Despite  the  knowledge  that  the  existence  of  a  pitching  mode 
would  affect  the  accuracy  of  the  data,  the  firing  was  conducted. 


FIG.  19.  Rocket  Motor  and  Stand  Displacement  as  a  Function 
of  Frequency  (Rocket  Motor  Carcass) . 
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LIVE  MOTOR  FIRING 


An  evaluation  test  of  the  experimental  system  with  the  NOTS- designed 
three-component  static  thrust  stand  was  conducted  with  the  firing  of  a 
large  rocket  motor  on  29  June  1962.  The  system  was  operated  as  a  straight 
sampling  type.  With  the  live  motor  mounted  in  the  stand  and  the  mass¬ 
measuring  system  in  operation,  a  natural  frequency  of  10.88  cps  was  indi¬ 
cated.  This  reflects  a  corrected  propellant  mass  (calculated  mass  minus 
the  mass  of  the  motor-stand  hardware)  of  approximately  466  slugs  or  a 
propellant  weight  of  15,000  lb.  At  the  end  of  burning  a  natural  fre¬ 
quency  of  16.88  cps  was  indicated- -which  reflects  a  motor-stand  hardware 
weight  of  approximately  10,000  lb. 

The  control  system  remained  locked- in  throughout  burning  and  the 
natural  frequency  of  the  spring-mass  system  during  burning  was  obtained. 
The  first  full-scale  demonstration  of  the  feasibility  of  this  technique 
for  determining  mass  during  burning  is  shown  in  a  plot  of  mass  versus 
burning  time  derived  from  the  natural  frequency  data  (Fig.  2l). 


TIME,  SEC 


FIG.  21.  Vibrating  Mass  as  a  Function  of  Burning  Time. 
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It  is  evident  from  Fig.  21  that,  except  for  the  period  from  5  to 
15  sec  during  "burning,  the  calculation  of  mass  using  natural -frequency 
information  is  generally  within  2$  of  the  mass  as  derived  from  internal 
ballistics  or  pressure  information.  The  graph  was  plotted  by  smoothing 
or  averaging  the  data  to  obtain  a  point  for  each  second  of  operation. 

Several  anomalies  are  apparent  in  the  5-  to  15-sec  period  of  burn¬ 
ing.  Analysis  of  the  data  of  the  front  and  rear  displacement  sensors 
indicates  that  the  pitching  motion  is  predominately  active  during  this 
state  of  operation  since  the  natural  frequency  of  the  pitching  mode  is 
equal  to  that  of  the  translational  mode  during  this  period.  It  is  ap¬ 
parent  that  the  control  system  is  attempting  to  lock-on  the  pitching 
mode  and  is  thus  switching  between  the  pitching  and  translational  motions. 
The  natural  frequency  of  the  pitching  mode  and  the  natural  frequency  of 
the  translational  mode  increase  with  a  decrease  in  mass  but  the  pitching 
mode  frequency  increases  at  a  much  lower  rate.  Therefore,  after  15  sec 
of  burning,  the  natural  frequency  of  the  pitching  mode  is  sufficiently 
separated  from  the  desired  translational  natural  frequency  so  that  its 
interference  with  the  control  system  is  minimized.  Mass  calculations 
from  this  point  to  the  end  of  burning,  using  frequency  data,  are  gener¬ 
ally  within  2 $  of  the  mass  plotted  from  the  internal  ballistics  data. 

The  pitching  motion  contributes  additional  error  since  the  magni¬ 
tudes  of  the  displacement  sensors,  which  are  located,  over  the  front  and 
rear  springs,  are  not  equal.  Since  the  waveform  representing  system  re¬ 
sponse  is  obtained  by  summing  the  sinusoidal  outputs  of  the  sensors,  it 
'has  a  phase  relationship  dependent  upon  displacement  magnitude.  The 
extent  to  which  this  adds  to  frequency  error  is  not  known  since  true- 
system  response  is  not  readily  available  from  existing  data;  however,  in 
a  spring-mass  system  with  a  very  small  damping  ratio  (approximately  0.1$), 
a  relatively  large  phase  error,  say  10°,  would  result  in  an  approximate 
1$  error  in  the  natural  frequency  and  an  approximate  2$  error  in  mass 
calculation. 


Another  factor  associated  with  the  pitching  motion  which  contrib¬ 
utes  directly  to  error  in  calculations  of  mass  is  the  assumption  that 


M 


Kf  +  Kr 


(51) 


This  equation  describes  a  spring-mass  system  with  a  single  degree  of 
freedom.  To  derive  the  approximate  equation  including  the  pitching  mode, 
consider  a  spring-mass  system  with  little  damping,  wherein: 


PE 

KE  ~ 


1 
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i(K.+K  )  z2  +  A(i?K_  +  i2K  )  e 
*  f  r  m  rf  r  r'  : 


|  M(o)  Z  )  +  i  j(co  0  )‘ 

2  v  n  m'  2  v  n  irr 


=  1 


Solving  for  mass: 


M  = 


(Kf  +  k  )  z2  +  (j!k-  +  A  )  e2  je2 


m 


f  f 

2  2 
CD  Z 
n  m 


r  r' 


m 


m 


m 


(52) 


where 


PE  =  maximum  potential  energy,  lb-ft 

KE  =  maximum  kinetic  energy,  lb- ft 

2 

J  =s  inertia,  lb-ft- sec 


£  =  distance  from  c,.g.  to  forward  support,  ft 
i  =  distance  from  c.g.  to  rear  support ,  ft 
K  =  front  spring  constant,  lb/ ft 
=  rear  spring  constant,  lb/ ft 

Cf  =  front  damping,  lb~^e- 

-  ,  .  lb-sec 

C  =  rear  damping,  — — — 
r  it 

Z  =  maximum  vertical  displacement,  ft 

0)  =  natural  frequency  vertical  plane,  rad/sec 

9^  =  maximum  pitching  angle,  rad 


A  rough  calculation  using  quick-look  analog  data  for  obtaining  0 
and  estimating  a  value  for  J  has  shown  that  an  error  of  2 $  in  mass  de¬ 
termination  is  probable  if  the  simplified  mass  equation  is  used  (see 
Fig.. 21  for  correction  of  the  two  end  points). 
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Another  possible  source  of 
spring  nonlinearity.  A  dynamic 


error  in  the  determination  of  mass  is 
spring  constant  can  be  derived  as  follows 


where 


“l 


mass  of  live  motor 


mass  of  motor  Carcass 


Rearranging,  the  equation  becomes 


M  K 

ml  =  — 

CD 


and 


n2 


and  combining  terms,  gives 


w  ,,  K  K 

-  “b  =  ~ — r 

CD  CD 

nl  n2 


Again  rearranging,  the  equation  becomes 


m 


) 


and  finally, 


K(ay) 


( 


m 


\  2  2 
x  CD  -  CD 

n 


n. 


1 


(53) 


Using  the  end  point  frequencies  and  the  change  in  mass  obtained  by 
weighing  the  rocket  motor  before  and  after  burning,  a  dynamic  spring  con¬ 
stant  of  3*72  x  106  lb/ ft  was  obtained.  This  figure  is  within  approxi- 
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mately  k . 5%  of  the  average  calculated  spring  constant  of  3*56  x  10°  lb/ft 
obtained  using  static  deflections  and  known  weights  in  the  mass  calcula¬ 
tions  shown  in  Fig.  21.  Although  the  pitching  mode  has  been  neglected 
in  calculating  for  this  magnitude  of  difference  suggests  that  a 

a  measure  of  nonlinearity  may  exist. 

Analysis  of  the  data  from  this  firing  in  correlation  with  the  analog 
computer  study  (Appendix  A)  indicates  that  the  discrepancy  in  the  calcula¬ 
tion  of  mass  from  actual  frequency  as  compared  with  internal  ballistics 
data  can  be  attributed  to  the  existence  of  the  pitching  mode  and  to  the 
momentary  loss  of  vertical  displacement. 

K 

It  is  estimated  that  using  the  simplified  equation,  mass  =  - 

oT 

n 

(which  neglects  the  pitching  mode)  contributed  an  error  of  approximately 
2 °!o  in  the  computation  of  mass.  It  should  be  noted,  however,  that  informa¬ 
tion  obtained  in  the  simulation  study  indicated  that  the  pitching  mode 
error  as  well  as  the  error  in  the  indicated  natural  frequency  may  be  due 
to  an  additional  form  of  damping  or  structural  phenomena. 

Experimental  tests  were  conducted  to  determine  what  types  of  damping 
existed  and  if  undesirable  structural  oscillations  were  present.  It  was 
found  that  structural  as  well  as  viscous  damping  existed  and  that  pitch¬ 
ing  was  being  introduced  by  coupling  through  damping  and  by  the  moments 
generated  by  the  structural  oscillation  of  the  thrust  take-out  member. 

Although  structural  damping  adds  to  the  nonlinearities  of  the  system, 
it  is  not  particularly  undesirable  in  system  operation  since,  in  the  trans 
lational  mode,  the  undamped  natural  frequency  is  independent  of  such  damp¬ 
ing.  It  was  observed  that  the  natural  frequency  of  the  thrust  take-out 
member  was  approximately  15  cps.  With  this  relatively  low  natural  fre¬ 
quency,  the  member  was  being  excited  through  the  oscillations  in  the  trans 
lational  mode.  A  solution  to  this  problem  might  be  to  stiffen  the  thrust 
take-out  member  (which  can  be  accomplished  quite  easily)  to  raise  its  nat¬ 
ural  frequency  and  thus  reduce  the  energy  transmitted  from  the  transla¬ 
tional  mode.  The  pitching  motion  further  introduced  error  in  the  phase 
relationship  of  the  displacement  waveform.  This  spurious  error  in  phase 
relationship  was  the  result  of  sensing  the  unequal  translational  motion 
of  the  rocket  motor  at  the  front  and  rear  supports  and  was  reflected  in 
the  control  system  as  true  freqency  error.  The  frequency  error  from  this 
source  (generally  spasmodic)  is  predominately  active  during  the  first  15 
sec  of  burning.  It  was  observed  that  the  pitching  motion  was  most  active 
during  this  period  when  the  natural  frequency  of  the  pitching  mode  was 
equal  to  and  traversed  by  the  continually  increasing  natural  frequency  of 
the  spring-mass  system. 

Momentary  loss  of  vertical  displacement  resulting  both  from  excessive 
phase  error  and  contributions  from  external  disturbances  caused  an  over¬ 
excitation  of  the  control  system.  In  this  state,  large  magnitudes  of  os¬ 
cillation  are  introduced  and  the  system  approaches  instability.  These 
oscillations  are  apparent  In  the  5"  to  15-sec  period  of  burning. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  theoretical  analysis  presented  in  this  report  gives  a  reasonable 
prediction  of  the  mass  measurement  system  performance.  As  well  as  provid¬ 
ing  component  design  criteria,  the  root-locus  method  becomes  a  convenient 
tool  for  field  optimization  of  the  system.  The  nature  of  the  root-locus 
method  is  such  that  the  time  constant  of  each  component  is  easily  obtained 
and  its  affect  on  the  system  can  be  determined.  The  affect  of  various 
types  of  compensating  circuits  (i.e.,  lead,  lag,  lead-lag,  or  lag-lead) 
on  the  root  locus,  and  hence  on  the  system,  can  also  be  conveniently  de¬ 
termined. 

The  assumption  that  the  spring-mass  nonlinear  transfer  function  can 
be  replaced  with  that  of  a  simple  RC  lag  network  is  not  theoretically  just¬ 
ified.  It  is  clear,  however,  from  the  analog  studies  that  a  lag  does  exist 
in  the  feedback  loop  due  to  sudden  changes  in  mass.  It  is  also  clear  that 
the  time  response  of  the  lag  changes  between  the  two  end  points  of  opera¬ 
tion.  This  is  probably  due  to  a  change  in  system  damping  between  these 
operating  points.  A  detailed  analysis  of  this  lag  is  described  in  NAVWEPS 
Report  77^1  (see  Bibliography).  This  report  describes  an  electrical  ana¬ 
log  of  the  spring-mass  system  which  is  designed  and  constructed.  The 
breadboard  model  is  laboratory  tested  to  determine  the  lag  characteristics. 

The  experimental  system,  through  the  successful  firing  of  a  large 
rocket  motor,  has  demonstrated  the  feasibility  of  this  technique  for  de¬ 
termining  the  mass  of  large  rocket  motors  during  burning. 

The  steady- state  performance  characteristics  indicate  that  the  system 
(straight  sampling  type)  is  capable  of  maintaining  a  translational  oscil¬ 
lation  within  1$  of  the  undamped  natural  frequency  of  the  spring-mass  sys¬ 
tem.  A  1  $  frequency  error  in  a  single  degree  of  freedom  system  will  cor¬ 
respond  to  an  approximate  2$  error  in  mass  measurement.  Since  the  present 
methods  of  mass  determination  are  considered  to  have  an  accuracy  of  5 1°> 
it  is  obvious  that  the  new  technique  will  improve  accuracies  significantly. 
The  main  limitation  of  the  experimental  system  is  its  inability  to  respond 
rapidly  to  sudden  changes  in  mass  burning  rates.  It  is  estimated  that 
sudden  mass  changes  in  the  order  of  1  $  of  the  total  vibrating  mass  could 
be  detected.  Detailed  information  on  the  mass  burning  rates  during  a  period 
of  irregular  burning  cannot  be  obtained  because  of  the  smoothing  or  inte¬ 
grating  characteristic  of  the  control  system. 

It  is  pointed  out  that  this  experimental  system  was  developed  primar¬ 
ily  to  demonstrate  the  feasibility  of  this  technique  of  mass  measurement. 
Considerable  improvements  can  be  incorporated  into  an  operational  system 
to  improve  steady-state  and  transient  characteristics.  Several  of  these 
improvements  are  discussed  in  the  following  paragraphs. 
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OPERATIONAL  SYSTEM 

As  noted  above,  the  design  of  an  operational  control  system  involves 
two  prime  considerations:  minimization  of  steady-state  error  and  maximi¬ 
zation  of  frequency  (transient)  response.  The  experience  gained  in  the 
successful  operation  of  the  experimental  system  has  provided  considerable 
insight  in  the  selection  of  components  to  adequately  meet  these  two  major 
objectives.  Following  are  several  suggested  improvements  for  each  of  the 
major  experimental  system  components. 


Force  Generator 


The  force  generator  used  in  the  experimental  system  was  a  mechanical 
device  with  rotating  parts  (Fig.  5)*  Although  the  device  generated  a 
fairly  smooth  sinusoidal  waveform  some  backlash  did  exist.  This  backlash 
appeared  minimal  and  was  not  considered  in  the  theoretical  anaylsis;  how¬ 
ever,  the  existence  of  backlash  is  certainly  undesirable.  The  force  gen¬ 
erator  was  driven  with  a  DC  motor  which  in  turn  was  controlled  by  a  mag¬ 
netic  amplifier.  This  relatively  long  power  train  results  in  excessive 
response  lags  and  ultimately  limits  the  frequency  response  of  the  closed- 
loop  system.  As  an  improvement  for  this  portion  of  the  system  it  is  rec¬ 
ommended  that  an  electrodynamic  shaker  be  employed.  Electrodynamic  shak¬ 
ers  that  will  fulfill  frequency  and  force  output  requirements  are  commer¬ 
cially  available.  In  addition  to  the  advantage  of  generating  a  smooth 
undistorted  waveform,  the  main  improvement  the  electrodynamic  shaker  pro¬ 
vides  is  a  large  reduction  in  response  time.  It  is  estimated  that  the 
response  times  associated  with  this  portion  of  the  experimental  system 
could  be  reduced  by  a  factor  of  10  with  the  selection  of  a  suitable  elec¬ 
trodynamic  shaker.  These  shakers  are  available  with  both  velocity  and 
acceleration  feedback  so  that  a  closed-loop  system  similar  to  that  of,  the 
experimental  system  could  be  used. 


Phase  Comparator 

The  phase  comparator  used  in  the  experimental  system  was  a  sampling 
device  in  that  the  phase  error  was  determined  at  each  zero  crossover  of 
the  displacement  wave.  It  is  apparent  that  an  improvement  in  system  oper¬ 
ation  could  be  realized  if  a  phase  comparator  were  available  for  a  con¬ 
tinuous  comparison  of  the  phases  of  two  sinusoidal  waveforms.  The 
complication  is  in  the  design  of  a  phase  comparator  capable  of  comparing 
the  phase  relationship  of  two  waveforms,  one  of  which  exhibits  a  variable 
amplitude.  Some  methods  of  low  frequency  phase  detection  are  described 
in  a  paper  by  T.  F.  Bogart,  Jr.* 


Technical  paper  dated  January  1963  written  to  fulfill  thesis  require¬ 
ment  for  a  Masters  Degree  at  the  University  of  California,  Los  Angeles. 


WAVWEFS  REPORT  835**- 


Another  method  of  improving  this  portion  of  the  control  system  would 
be  to  use  a  sampling  phase  comparator  as  was  used  in  the  experimental  sys¬ 
tem  and  increase  the  sampling  rate.  Although  the  following  methods  have 
not  been  investigated  in  detail,  this  effect  might  be  accomplished  through 
the  use  of  frequency  multipliers  of  electronic  type  or,  in  the  case  of 
the  sine-cosine  potentiometer,  through  a  simple  gearing  arrangement. 


Spring-Mass  Dynamics 

Although  this  area  has  been  touched  only  lightly  in  this  report  it 
remains  a  very  important  element  in  the  successful  use  of  this  technique 
of  mass  measurement.  The  static  thrust  stand  must  be  designed  to  meet 
special  requirements.  The  stand  must  be  adjusted  in  the  field  to  elimin¬ 
ate  entirely  if  possible,  or  to  at  least  minimize,  pitching.  This  would 
enable  the  use  of  single  degree  of  freedom  theory  and  much  simpler  oper¬ 
ation.  NAVWEPS  Report  8353*  describes  the  design  criteria  for  construc¬ 
tion  of  thrust  stands  to  be  used  with  this  technique  of  mass  measurement. 


1  jo  System 

It  is  felt  that  modification  of  the  control  system  described  in  this 
report  to  incorporate  the  use  of  an  electrodynamic  shaker  and  by  careful 
thrust- stand  design,  a  system  can  be  developed  which  will  provide  for  the 
determination  of  mass  to  an  accuracy  of  better  than  1  $  and  will  enable 
the  extraction  of  data  on  instantaneous  burning  rates  during  periods  of 
irregular  burning. 


U.  S.  Naval  Ordnance  Test  Station.  Design  Criteria  for  Large 
Accurate  Solid-Propellant  Static-Thrust  Stands,  by  D.  P.  Ankeney  and  C. 
E.  Woods.  China  Lake,  Calif.,  NOTS,  June  1963*  (NAVVEPS  Report  8353; 
WOTS  TP  32^0). 
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Appendix  A 


ANALOG  SIMULATION 


Data  analysis  of  the  first  evaluation  firing  using  the  mass  measur¬ 
ing  system  pointed  out  some  apparent  anomalies  in  the  data  and  in  the 
system's  overall  performance.  In  order  to  better  define  system  behavior 
and  to  predict  minor  in-house  modifications  which  could  be  made  to  the 
experimental  system  before  conducting  the  second  evaluation  firing*  an 
analog  computer  simulation  of  the  complete  system  was  made.  In  addition, 
an  investigation  to  determine  system  performance  under  jetavator  reaction 
was  conducted. 


ANALOG  SETUP 


The  differential  equations  describing  the.  spring-mass-force  system 
were  written  in  two  degrees  of  freedom  (Appendix  B)  taking  the  vertical 
motion  of  the  c.g.  of  the  rocket  motor  and  the  pitching  motion  of  the 
rocket  motor  about  its  c.g.  into  consideration.  The  equations  were  ex¬ 
tended  to  include  both  front  and  rear  vertical-motion  sensors  (potentiom¬ 
eters).  To  simulate  the  two-degree  of  freedom  experimental  system,  the 
signals  from  the  sensors  were  summed  to  provide  one  composite  signal  rep¬ 
resenting  the  vertical  displacement  of  the  rocket  motor.  A  schematic  of 
the  analog  simulation  is  presented  in  Appendix  B. 

In  setting  up  the  analog  computer  study  it  was  felt  that  consider¬ 
able  valuable  information  (in  addition  to  that  defining  general  system 
characteristics)  would  be  obtained  if  the  experimental  system  were  very 
accurately  simulated.  For  instance,  the  same  improved  system  performance 
produced  in  the  analog  simulation  by  small  changes  in  electrical  charac¬ 
teristics  and  minor  hardware  modifications  could  also  be  expected  from 
the  experimental  system  by  making  the  same  component  changes. 

The  procedure  used  in  attempting  to  set  up  the  analog  computer  to 
simulate  the  experimental  system  accurately  was  to  adjust  the  basic 
(experimental  or  calculated)  values  of  various  coefficients  in  the  dif¬ 
ferential  equations  that  were  known  or  believed  to  change  during  burning 


The  second  firing  was  not  conducted  because  of  a  rocket-motor  mal¬ 
function  which  resulted  in  partial  destruction  of  the  drive  motor  and 
shaker . 
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and  then  try  to  retrace  the  vertical  displacement- curve  data  shown  in 
Figs.  19  and  20.  To  obtain  data  points  for  the  curves,  the  system  was 
operated  open  loop  and  the  front  and  rear  rocket-motor  displacement  was 
observed  as  the  oscillation  frequency  was  varied.  The  displacement 
curves  of  the  final  analog  setup  are  shown  in  Figs.  22  and  23*  Figures 
19  and  22  are  curves  for  a  rocket-motor  carcass  mounted  in  the  stand  and 
Figs.  20  and  23  are  curves  for  a  live  rocket  motor  mounted  in  the  stand. 

In  order  to  reproduce  the  curves  plotted  in  Figs.  19  and  20  for  the  ana¬ 
log  simulation,  it  was  necessary  to  use  time-varying  functions  for  rocket- 
motor  inertia,  system  damping,  and  all  terms  associated  with  a  shift  in 
rocket-motor  c.g. 


FIG.  22.  Rocket  Motor  and  Stand  Displacement  as  a  Function 
of  Frequency  (Rocket  Motor  Carcass),  Analog  Simulation. 


kh 
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FIG.  23.  Rocket  Motor  and  Stand  Displacement 
as  a  Function  of  Frequency  (Live  Rocket  Motor) 
Analog  Simulation. 


The  controlled  functions  were  made  to  vary  linearly  with  time  from 
rocket-motor  ignition  to  "burnout.  Curves  of  each  function  versus  time 
appear  in  Appendix  B.  Although  the  displacement  curves  are  only  fairly 
comparable  to  those  of  the  experimental  system,  they  represent  the  best 
simulation  obtainable  with  the  present  differential  equations.  Addi¬ 
tional  terms  would  be  required  to  improve  the  accuracy  of  simulation. 
Structural  damping,  which  is  amplitude  dependent,  has  been  suggested 
as  a  highly  probable  required  term- -structural  damping  was  detected  in 
later  experimental  tests. 

The  actual  differential  equations  solved  include  viscous  damping 
only.  In  the  final  analog  setup  it  was  necessary  to  double  the  system 
damping  from  ignition  to  burnout.  It  was  further  necessary  to  reverse 
the  distribution  of  the  system  damping  between  front  and  rear  so  that 
the  damping  in  the  rear  exceeded  that  in  the  front.  This  was  necessary 
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in  order  to  reproduce  the  pitching  evident  in  the  experimental  data. 

It  has  since  been  determined  that  the  thrust  take-out  member  was  also 
contributing  to  the  pitching  motion. 

Since  the  differential  equations  solved  include  time-varying  coeffi¬ 
cient  s,  they  are  considered  nonlinear  and  their  degree  of  similarity  to 
the  experimental  system,  as  compared  by  vertical  displacement  curves,  is 
difficult  to  ascertain.  It  was  felt  however,  that  the  major  character¬ 
istics  of  the  displacement  curves  were  reproduced  and  that  reasonable 
similarity  existed. 

To  further  compare  the  analog  simulation  and  the  experimental  sys¬ 
tem,  all  of  the  components  associated  with  the  experimental  system  were 
simulated  and  tested.  Response  curves  were  obtained  by  monitoring  the 
frequency  of  the  driving  function  (motor  speed  0)  and  closing  the  con¬ 
trol  loop  while  the  system  was  oscillating  out  of  resonance.  The  re¬ 
sulting  curve  indicates  the  system  response  to  a  step  input  in  phase 
error  (an  approximated  step  change  in  mass).  Phase  error  is  defined  as 
the  phase  difference  in  degrees  between  the  forcing  function  and  the 
rocket-motor- stand  response  function  minus  90° ,  or: 


where  Qe  is  the  instantaneous  phase  error,  Of  is  the  phase  angle  of  the 
forcing  function  with  respect  to  )me  reference,  and  0r  is  the  phase 
angle  of  the  rocket-motor-stand  response  function  with  respect  to  the 
same  reference.  The  (-90°)  is  required  since  Of  is  leading  0r  by  9®°  at 
the  natural  frequency.  The  analog  simulation  response  curves  for  the 
step  input  in  phase  error  with  a  rocket -motor  carcass  and  with  a  live 
rocket  motor  mounted  in  the  stand  are  shown  in  Figs.  24  and  25*  Figures 
17  and  18  are  curves  from  the  same  test  conducted  experimentally  with 
the  prototype  system.  The  forward-loop  gain,  the  feedback-loop  gain,  and 
the  compensating  networks  of  the  control  system  were  adjusted  in  an  ef¬ 
fort  to  improve  the  similarity  between  the  two  sets  of  curves.  Figures 
24  and  25  represent  the  best  simulation  attainable  using  constant  con¬ 
trol  system  characteristics;  i.e.,  closed-loop  gain  and  compensating  net¬ 
works  . 

It  was  noted  that  much  better  similarity  could  be  obtained  (when 
either  the  rocket-motor  carcass  or  the  live  rocket  motor  was  mounted 
in  the  stand)  if  the  control  system  characteristics  were  adjusted  under 
each  condition.  This  adjustment,  however,  cannot  presently  be  made  in 
the  field,  and  it  was  decided  to  continue  the  simulation  study  using  the 
analog  setup  indicated  above  rather  than  to  attempt  to  incorporate  time- 
varying  functions  in  the  control  system. 

Considerable  difference  is  apparent  between  the  two  sets  of  response 
curves.  The  degree  of  dissimilarity  between  the  analog  setup  and  the  ex¬ 
perimental  system  is  such  that  the  significance  of  the  analog  study  will 
be  limited  to  applications  of  a  general  nature. 
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FIG.  24.  System  Response  as  a  Function  of  Time  to  a  Simu¬ 
lated.  Step  Input  (Rocket  Motor  Carcass),  Analog  Simulation. 


FIG.  25.  System  Response  as  a  Function  of  Time  to  a  Simu¬ 
lated  Step  Input  (Live  Rocket  Motor),  Analog  Simulation. 
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COMPUTATION  OF  MASS 


The  curve  of  Fig.  2 6  is  mass  versus  time  and  represents  the  simu¬ 
lated  linearly  changing  mass.  The  curves  of  Fig.  27  are  mass  versus  time 
and  were  obtained  under  closed-loop  operation  by  computing  mass  from  con¬ 
tinuous  knowledge  of  the  driving- function  frequency  and  the  spring  con¬ 
stant  of  the  spring-mass  system. 

It  is  noted  that  very  little  difference  exists  between  the  computed 
curves  of  mass  versus  time  with  and  without  the  pitching  mode  terms. 

That  is,  the  computed  curves  of  mass  versus  time  obtained  when  simulat¬ 
ing  a  two- degree -of- freedom  spring-mass  system  are  essentially  retraced 
when  the  spring-mass  system  is  simulated  in  one  degree  of  freedom.  The 
one-degree-of- freedom  simulation  was  accomplished  by  eliminating  all 
terms  in  Eqs.  (49)  and  (50)  associated  with  the  pitching  angle  f*  The 
resulting  equation  is 


M(t)z  +(0^  +  0  )(t)Z  +  (X  +  K  )z  =  m.to2sina)  t 
\  /  n  f  r  '  v  f  r  o  o 

An  approximate  2$  error  in  the  mass  computation  due  to  rocket  motor 
pitching  was  expected  since  this  was  the  error  estimated  in  the  live  mo¬ 
tor  test.  This  error  is  not  apparent  in  the  analog  simulation.  An  ap¬ 
preciable  amount  of  pitching  was  evident  in  the  analog  setup  which  tends 
to  indicate  that  pitching  does  not  contribute  greatly  to  frequency  error. 

It  should  also  be  pointed  out  that  it  was  not  possible  to  excite 
the  pitching  mode  when  simulating  a  statically  decoupled  system.  That 
is,  when  setting  the  term  (ifKf  -  irKr)  =  0.  This  tends  to  indicate  that 
damping  other  than  viscous  is  present  or  perhaps  that  coupling  is  exis¬ 
tent  through  a  structural  difficiency  or  through  the  shift  in  c.g.  It 
was  necessary  to  introduce  a  relatively  large  coupling  term  (i^Kf  -  irKp) 
in  order  to  cause  pitching  and  enable  the  retracing  of  the  displacement 
curves  of  Figs.  22  and  23 •  The  overall  accuracy  in  analog  mass  computa¬ 
tion  is  improved  over  the  indicated  results  of  the  experimental  system 
in  the  evaluation  firing.  The  error  in  mass  computation  is  approximately 
2 *fo  at  ignition  and  decays  to  essentially  zero  by  the  end  of  burning. 
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ulation. 


FIG.  27*  Vibrating  Mass  as  a  Function  of  Burning  Time 
Computed  From  Frequency  Data,  Analog  Simulation. 
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JETAVATOR  REACTION 


The  curves  of  Figs.  28  and  29  are  computed  mass  versus  time  with 
the  system  "being  subjected  to  jetavator  reaction  of  various  frequencies 
and  waveforms  (Appendix  B).  It  was  found  that  the  system  as  simulated 
was  essentially  unaffected  (actually  performance  was  improved)  by  a  sinus¬ 
oidal  jetavator  reaction  frequency  of  1  or  2  cps.  The  system  became  un¬ 
stable  as  the  jetavator  reaction  frequency  approached  10  cps  (sinusoidal). 
In  each  case  where  instability  was  observed,  it  was  noted  that  the  con¬ 
trol  loop  remained  unaffected  and  that  instability  took  place  through 
the  progressive  buildup  of  the  vertical  displacement.  The  system  again 
became  unstable  with  the  application  of  a  square  wave  jetavator  reaction 
of  1-cps  frequency.  The  system  was  essentially  unaffected,  however,  when 
subjected  to  a  triangular  waveform  with  a  frequency  of  2  cps. 


ROCKET-MOTOR  LAG 


From  data  obtained  with  the  experimental  system,  it  is  apparent 
that  a  lag  relative  to  changes  in  driving- function  frequency  exists  in 
the  spring-mass  system.  It  was  not  known  if  this  lagging  effect  would  be 
sufficient  to  hinder  closed-loop  performance  under  conditions  of  normal 
or  irregular  propellant  burning.  An  investigation  was  conducted  on  the 
analog  system  to  better  define  the  rocket-motor  lag  characteristics. 

The  investigation  consisted  of  effecting  step  changes  in  mass  while  hold¬ 
ing  the  oscillation  frequency  of  the  driving  function  constant  and  ob¬ 
serving  the  phase  angle  between  the  spring-mass  system  and  the  driving 
function.  The  curves  in  Fig.  30  show  phase  angle  versus  time  for  the 
two  conditions;  i.e.,  live  and  carcass  rocket  motors  mounted  in  the  stand. 
Comparison  of  the  response  curves  of  Figs.  24  and  25  with  those  of  Fig.  30 
shows  that  the  response  times  are  approximately  the  same  in  each  case. 
While  this  indicates  that  the  rocket-motor  lag  in  a  spring-mass  system 
might  very  well  affect  overall  system  performance,  it  is  noted  that  the 
simulated  system  did  perform  satisfactorily  with  the  lag. 

Although  the  experimental  system  is  only  generally  simulated  and  the 
rocket-motor  lag  curves  do  not  necessarily  reflect  field  conditions,  it 
is  reasonable  to  assume  that  this  lag  does  not  seriously  affect  the  per¬ 
formance  of  the  present  experimental  system.  However,  if  a  degree  of 
accuracy  better  than  the  presently  feasible  2$  is  to  be  attained,  addi¬ 
tional  investigation  in  this  area  may  be  necessary. 
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A  ROCKET  MOTOR  CARCASS  PHASE  DIFF  IS  THE  ANGLE  AT  WHICH  THE 

(DATA  GOOD  TO  GO*  ONLY)  MASS  SYSTEM  LAGS  THE  FORCING  FUNCTION. 

_  i  ii/r  unrfto  THIS  PLOT  WAS  OBTAINED  BY  INITIATING  A 

o  LIVE  MOTOR 

STEP  CHANGE  IN  MASS  AT  TIME'O  ANO 
HOLDING  THE  FREQUENCY  OF  OSCILLATION 
CONSTANT  (OPEN  LOOP) 

STEP  CHANGE'  40  SLUGS 


TIME,  SEC 


FIG.  30.  Phase  Angle  Response  of  Spring-Mass  System  as  a 
Function  of  Time,  Analog  Simulation. 


SUMMARY 


Analog  computer  simulation  of  the  experimental  mass  measurement  sys¬ 
tem  was  undertaken  with  the  idea  that  the  system  could  be  accurately  sim¬ 
ulated  and  that  reproduction  of  experimental  data  could  be  made.  Detailed 
changes  could  then  be  made  on  the  analog  setup  to  study  and  improve  system 
performance  and,  by  making  the  appropriate  changes  on  the  experimental 
system,  the  same  results  could  be  expected. 

The  effective  similarity  between  the  final  analog  setup  and  the  ex¬ 
perimental  system  was  found,  however,  to  be  questionable  since  only  part 
of  the  data  obtained  with  the  experimental  system  could  be  reproduced  in 
the  simulation.  The  data  reproduced  are  of  vertical  displacement  versus 
frequency  for  the  two  conditions  of  a  rocket-motor  carcass  and  a  live 
rocket  motor  mounted  in  the  thrust  stand.  In  order  to  obtain  correlation 
with  experimental  data  it  was  necessary  to  vary,  linearly  with  time  dur¬ 
ing  burning,  several  of  the  coefficients  of  the  differential  equations 
solved.  The  extent  to  which  some  of  these  coefficients  had  to  be  varied 
does  not  appear  to  be  practical.  The  conclusion  is,  then,  that  the  analog 
setup  must  be  considered  in  general  terms  only  and  that  it  cannot  be  used 
for  detailed  design  studies.  Since  the  analog  setup  is  of  the  same  type 
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of  system  as  the  experimental  system  (i.e.,  velocity- controlled,  phase- 
sensitive)  with  identical  type  components,  it  is  reasonable  to  speculate 
that  the  general  system  characteristics  are  also  similar. 

It  should  be  noted  that  several  tests  conducted  with  the  analog 
setup  have  indicated  that  damping  is  a  critical  factor  in  depicting  the 
similarity  between  the  analog  setup  and  the  experimental  system.  These 
tests  support  a  previous  theory  that  structural  (amplitude -dependent)  as 
well  as  viscous  (velocity- sensitive)  damping  is  existent  in  the  experi¬ 
mental  system.  It  is  reasonable  to  assume  that  improved  correlation  could 
be  obtained  in  a  future  study  by  using  additional  damping  terms.  Appro¬ 
priate  terms  would  also  be  required  to  simulate  the  pitching  moment  con¬ 
tributed  by  the  thrust  take-out  member. 

The  following  conclusions  are  based  on  the  preceding  discussion  and 
are  applied  to  the  experimental  system  in  terms  of  gross  information  only. 

The  vertical  frequency  characteristic  of  the  closed-loop  system 
(analog  setup)  to  a  linearly  changing  mass  is  one  of  controlled  symmetri¬ 
cal  oscillation  about  the  mean  or  natural  frequency.  The  computer  set¬ 
up  did  not  produce  excessive  oscillation  during  the  5-  to  15-sec  period 
of  burning  as  did  the  experimental  system.  This  is  attributed  to  the 
need  for  additional  terms  in  the  differential  equations  to  cause  the  ex¬ 
cessive  pitching  apparent  in  the  evaluation  test  during  this  period  of 
burning. 

A  rocket-motor  lag  to  a  step  change  in  mass  was  apparent  in  the  spring 
mass  system.  The  effect  of  such  a  lag  on  the  system’s  performance  is  small 
however,  and  should  not  prevent  mass  determination  to  an  accuracy  of  2 °jo* 

Jetavator  reaction  frequencies  of  2  cps  sinusoidal  and  triangular 
did  not  adversely  affect  system  performance.  The  system  was  unstable, 
however,  when  subjected  to  a  jetavator  reaction  sinusoidal  frequency  of 
10  cps  and  a  square  wave  frequency  of  1  cps. 

The  limited  computer  time  available  did  not-  permit  an  investigation  of 
system  component  characteristics  in  relation  to  their  affect  on  system  per¬ 
formance  nor  an  investigation  of  system  performance  when  operated  as  a  re¬ 
lay-sampling  system  (addition  of  relay  in  forward  loop).  Detailed  informa¬ 
tion  of  such  tests  would  not  be  direcly  applicable  to  the  experimental 
system  since  the  analog  setup  did  not  exactly  simulate  the  experimental 
system. 
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Appendix  B 

DERIVATION  OF  EQUATIONS  FOR  ANALOG  SIMULATION 
AND  ANALOG  SCHEMATIC  DIAGRAM* 

SPRING -MASS  SYSTEM  EQUATIONS 


F 10,  000  sin  (jO\t  FQ  =  mQlQ  t0o2  sin<0ot 


—  /)  -  ^ 

"  1 m  A  r  ~ 

- X  f 

Jcr  Cf  l_j 

j 

7’f7y  //  /  // // // / 


Rocket  Motor  Mounted  in  Stand 
Assume  two  degrees  of  freedom 

a.  Vertical  motion  of  c.  g. 

b.  Pitching  motion  about  c.  g. 

FQ  =  Force  generated  by  mechanical  shaker 
F ]_  =  Force  generated  by  jetavator 


From  F  =  ma 

MZ  =  [  m£  u>0s  sin  0)ot  -  (Kf  +  Kr)  Z  -  (Cf+Cr)  Z  -  (l  rKr-  t£Kf)  lj)  - 

(5^) 

(-t-rCr  -  -tfCf)  )p  -  10,  000  sin  0)1*] 

From  T  =  Ja. 

H'=  [ “('t-f 2  Kf  +  lraKr)4>-  (l£2C{+ lr3 CT)  il)-(lTKT  -If  Kf)  Z  - 

(55) 

(£rCr  ”^f  Cf)  Z  -  (lr  -f  1)  10,  000  sin  to  l*  ] 

- ^  — -  - — 

For  symbols  and  variable  coefficients  used  in  this  appendix,  see 
pages  57  and  59  j  the  analog  schematic  diagram  appears  on  page  60. 
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MECHANICAL  OSCILLATOR  EQUATION* 

Derivation  of  the  equation  of  torque  reflected  into  drive  motor 
mounted  on  inertial  axis. 


Unit  Vectors 

ey>  eZ>  ex 

then: 

&  =  e  z  £  sin  8  +  l  cos  d 

(56) 

and 

— J 1  d  — L  d 

t  -  ez  &  sin  6  +  e  y  cos  8 

(57) 

also 

Z  =  eZ  (  Z  t  Sin  0) 

(58) 

and 

*  =  °zlt  (z+<-sin  9) 

(59) 

—k  Lk  _ 1. 

r  =  Z  +  £ 

(60) 

"r*  =  "Z*  +  T 

(61) 

Torque  =  ^  (angular  momentum) 

L  X  dt 

(62) 

*This 

material  was  reproduced  from  notes  compiled  by  D. 

P.  Ankeney. 
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h  =  r  Xm  v 

(63) 

but 

dr^  _k 

dt  =  v  =  r 

(64) 

so: 

h  =  r  X  m  r 

(65) 

sub:  60-61  in  65 

h  =  (Z  +t)  Xm(Z  +t) 

(66) 

sub:  56-57-58-59 

in  66. 

h  —  m  ^0-y'C  cos  04* 

ez  (Z+^sin  0)]  X  cos0+^(Z+  l  sin  6)| 

(67) 

Expressing  6 2  in 

determinate  form 

ex  ®y  &Z 

L  f  =  mTt 

O  l  cos  0  Z  +  l  sin  0 

O  t|cos  S  A  (Z  +  t  Sin  0) 

(68) 

Solution  of  determinate  yields 

T'“i[ 

t,COS  6  ^  (Z  +  ^.sin0  )  -  t-^cos  0(Z  +  <£,sin0)j 

(69) 

or 

T  = 

m  -t  £  Z  0  sin  0  +  •{,  0  +  (Z  +  Z08)  cos  fij 

(70) 

MOTOR  EQUATIONS 

Separately  excited  shunt  motor 

Em  =IaRa  +  *1  6 

(71) 

T  =  Kq  Ia  (developed  torque) 

(72) 

T  =  J04B0  +  T  £  (required  torque) 

(73) 
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JETAVATOR  EQUATION 

F i  -  10,000  sin  oqt  (7^) 

Where  is  varied  with  following  waveforms: 

sinusoidal  1  to  10  cps 
triangular  1  and  2  cps 
square  1  cps 

DEFINITION  OF  SYMBOLS 

B  -  Viscous  damping  coefficient  (motor) 

Or  &  Of  =  Viscous  damping  coefficient  ( spring-mas  s  -  system) 

Em  =  Terminal  voltage  (drive  motor) 
la  -  Armature  current  (drive  motor) 

J  =  Polar  moment  of  inertia  (rocket  motor) 

J  =  Effective  moment  of  inertia  (drive  motor) 

Kl  =  Motor  constant  (drive  motor) 

Kr  &  Kf  =  Spring  constant  ( spring-mas s -system) 
l  =  Effective  moment  arm  (mechanical  oscillator) 

m  =  Effective  mass  (mechanical  oscillator) 

M  =  Sprung  weight  ( spring-mas s -system) 

Ra  =  Armature  resistance  (drive  motor) 

=  Angle  of  rotation  (pitching  mode) 

=  Angular  velocity  (pitching  mode) 

=  Angular  acceleration  (pitching  mode) 


* 

* 

•  » 
4> 
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&  =  Vertical  displacement  (rocket  motor) 

2  =  Vertical  velocity  (rocket  motor) 

£  =  Vertical  acceleration  (rocket  motor) 

8  =  Angle  of  rotation  (drive  motor) 

6  =  Angular  velocity  (drive  motor) 

6  =  Angular  acceleration  (drive  motor) 

h  =  Angular  momentum 
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time  scaled  ev  10 


m  M2.[-(Kf  *K,)I  -<C,«C,)Z  -U,  Kr  -  If  K|  )f -(JrCr.-itd()f 
♦"»*“*  SIN  Ml -10,000  SIN  w,t  { J€  T AVAT OftS  )J 

-<Jf*  1)10,000  SIN»(  HJCTAMATOAtjj 

(3)  Tf  *m/[i8*(i  +  2dI)C0S  •  +* SsiNtJ 

(4)  Erti*I„R0  *K,  9 

(5)  T*K|  l( 

(6) 

(7)  r«  10,000  SINWjl  ( JETAVATORS) 


amlog  schematic  diagram 
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Appendix  C 

DERIVATION  OF  SYSTEM  TRANSFER  FUNCTIONS 
AND  PHASE  COMPARATOR  SCHEMATIC 


INTEGRATOR 

A  high-gain  operational,  amplifier  was  used  in  the  following  config¬ 
uration  to  provide  the  integration  indicated  in  Fig.  7. 


The  transfer  function  for  an  operational  amplifier*  is  given  as 


_2  =  Z  f.b. 
e^  Z  in 

where 

Z  f.b.  =  feedback  impedance 
Z  in  =  input  impedance 
and  in  Laplace  notation 


z  f-b’  =  sc 

Z  in  =  2M 


1  x  10 


o 


then 


2S 


and 


K1  -5 


* 


Walter  J.  Karplus  and  Walter  W.  Sorous.  Analog  Methods .  New  York, 
McGraw-Hill,  1959-  P.  39- 

6l 


LEAD  COMPENSATION 

The  lead  compensation  indicated  in  the  forward  loop  of  Fig.  7  has 
the  following  configuration: 


Using  nodal  analysis: 


and 


*1  +  i2  +  *3  =  ° 


e 

o 

x3  =  1  M 


combining 


rearranging 


e 

o 


* 
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and 

eo  _  (S  +  l) 
e2  (S  +  11) 


in  the  same  manner  and  is 

f3  _  (S  +  1) 
eo  (S  +  50) 

from  which 

!i  _  (s  +  i)2 

e2  (S  +  11)  (S  +  50) 


—  is  obtained 
e 

o 


MINOR  LOOP  FEEDBACK 

The  differentiator  and  compensation  circuitry  indicated  in  the 
minor  loop  feedback  path  of  Fig.  7  utilizes  an  operational  amplifier 
connected  as  follows: 


-001  /if 


10  M 


— WV-  •  1 

IM 


>5.6K 


MAG -AMP 
BIAS  WINDING 
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Using  the  transfer  function  for  operational  amplifiers 

eo  =  Z  f «b. 
e  in  Z  in 

It  can  be  shown  that  the  feedback  capacitor  in  the  above  sketch  used  for 
filtering  the  differentiated  ripple  of  the  tachometer  output,  will  con¬ 
tribute  a  pole  in  the  overall  transfer  function  at  S  =  -100  and  as  such 
will  not  significantly  affect  the  root  locus  plot.  Therefore,  the  capac 
itor  will  not  be  included  in  the  following  derivation: 

el  =  10  M  =  20  S 

e  in  =  .5M  +  §  =  (_S  +  2) 

O 

and 

eo  10  M  10  (S  +  1) 

'i =  ~7  ‘  + 


s 


so  that  the  final  transfer  function  is 


e  in  e  in  e^ 


eo  _  200  S(S  +  1) 

6  in  =  (S  +  2 )2 

The  tachometer  constant  of  0.3  and  the  conversion  for  the  5*6K  load 
resistor  gives  a  final  transfer  function  of 

x7  10.8  x  10"3  S(S  +  1) 

9  (S  +  2)2 

K,  =  10.8  x  10 “3 
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PREAMPLIFIER 

The  preamplifier  indicated  in  Fig.  7  has  the  following  experimen¬ 
tally  determined  response  curve: 


U 


The  response  characteristic  is  considered  linear  over  the  necessary 
range  of  operation.  The  preamp  constant  is  found  to  be: 

^4  -3 

— ~  »  A,  =3x1 0  ^ 

1 


SIGNAL  SENSORS 

The  input  signals  to  the  phase  comparator  were  obtained  with  the 
use  of  carbon  film  potentiometers  (see  Fig.  31  for  approximate  hardware 
location) .  The  sinusoidal  signals  depicting  the  'driving -function  re¬ 
sponse  and  the  reference  response  were  obtained  by  attaching  a  sine- 
cosine,  continuous -turn  potentiometer  to  the  rotating  shaft  of  the  force 
generator.  The  signal  depicting  the  displacement  motion  was  obtained  by 
attaching  the  sliding  arm  of  a  linear -motion  potentiometer  to  the  thrust 
stand  leg.  Although  the  sine-cosine  unit  was  adequate  in  the  operation 
of  the  experimental  system,  these  units  have  a  relatively  short  life  and 
should  be  replaced  with  electromagnetic  devices  in  the  development  of  an 
operational  system. 
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FIG.  31*  Large  Rocket  Motor  in  Three - Component  Stand  Prior 
to  Firing. 


k 
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DC  DRIVE  MOTOR 

The  motor  constants,  Ra,  La*  and  K,  were  determined  experimentally 
using  the  procedure  described  in  the  literature  cited  below.*  Hie  values 
were  found  to  be: 


R  =  .86  ohms 
a 

L  =  4  n.h  (negligible) 
a 

T,  -  rc>  ohm-second 
K  —  1 « oo  « 

radian 

The  inertia  was  listed  by  the  manufacturer  for  the  motor  and  was  calcu¬ 
lated  for  the  force  generator.  The  total  inertia  is 

J  =  .71  ft-lb-sec2 

The  motor  time  constant  is  defined  as 


and  is  t  =  .1  seconds 
m 

Since  1  /%  positions  the  pole  in  the  motor  transfer  function,  the  follow¬ 
ing  expression  is  obtained,  neglecting  1^: 


iL  io 

eg  ~  S  (s  +  10) 


MAGNETIC  AMPLIFIER 

The  steady-state  response  curve  of  the  magnetic  amplifier  and  drive 
motor  was  determined  as  a  single  unit  under  operating  conditions,  see 
page  20.  The  transient  response  characteristic  of  the  magnetic  amplifier 
was  determined  experimentally  as  shown  in  the  following  sketch; 


Office  of  Chief  of  Ordnance,  Department  of  the  Army .  Servomechan¬ 
isms,  Section  k}  Power  Elements  and  System  Design.  August  1959*  (ORDP 

20-139). 
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OSCILLOSCOPE 

CAMERA 

RECORDING 


In  the  above  setup  the  fixed  resistance  Rf  is  connected  to  the 
field  output,  terminals  of  the  amplifier  to  simulate  the  field  of  the 
drive  motor.  The  resistance  simulates  the  motor  armature  load.  An 
appropriate  resistance  was  chosen  to  correspond  to  l/2  full  load.  With 
the  mag-amp  energized  the  switch  was  closed.  The  load  of  R2  is  then 
added  to  that  of  R^.  The  resulting  transient  was  recorded  and  is  inter¬ 
polated  below. 


STEP 

INPUT 

(VOLTAGE 

ACTUALLY 

DECREASES) 


SECONDS 


The  response  curve  approaches  that  of  an  exponential  resulting  from 
a  pure  R,  L  circuit  having  a  time  constant  of  0.1  seconds.  Thus,  the 
transfer  function  is 

e6  A2  _  A2 
*5  '  (s“+  Iprj  -  (s  +  10) 

is  determined  from  steady-state  characteristics  as  described  on 

page  20. 
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